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ABSTRACT 

Charmonium in Hot Medium. (December 2010) 

Xingbo Zhao, B.S., University of Science and Technology of China 

Chair of Advisory Committee: Dr. Ralf Rapp 

We investigate charmonium production in the hot medium created by heavy-ion 
colhsions by setting up a framework in which in-medium charmonium properties are 
constrained by thermal lattice QCD (IQCD) and subsequently implemented into ki- 
netic approaches. A Boltzmann transport equation is employed to describe the time 
evolution of the charmonium phase space distribution with the loss and gain term 
accounting for charmonium dissociation and regeneration (from charm quarks), re- 
spectively. The momentum dependence of the charmonium dissociation rate is worked 
out. The dominant process for in- medium charmonium regeneration is found to be 
a 3-to-2 process. Its corresponding regeneration rates from different input charm- 
quark momentum spectra are evaluated. Experimental data on J /ip production at 
CERN-SPS and BNL-RHIC are compared with our numerical results in terms of both 
rapidity-dependent inclusive yields and transverse momentum {pt) spectra. Within 
current uncertainties from (interpreting) IQCD data and from input charm-quark 
spectra the centrality dependence of J /ip production at SPS and RHIC (for both mid- 
and forward rapidity) is reasonably well reproduced. The J /ip pt data are shown to 
have a discriminating power for in-medium charmonium properties as inferred from 
different interpretations of IQCD results. 
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CHAPTER I 



INTRODUCTIOI 

The standard model has been estabhshed as the underlying theory of modern particle 
physics. It describes the three fundamental interactions (electromagnetic, weak and 
strong) in a unified framework. The standard model includes 12 spin- 1/2 particles 
(fermions). They are six quarks: up (m), down (d); charm (c), strange (s); top (t), 
bottom (6) and six leptons: electron (e), electron neutrino (ue), muon (/i), muon 
neutrino (i/^), tau (r), tau neutrino (z/,-). It also includes spin-1 particles (bosons) 
mediating the interaction between fermions. They are the photon mediating the elec- 
tromagnetic interaction, W^, W~ and Z bosons mediating the weak interactions and 
8 gluons mediating the strong interaction. The elementary particles in the standard 
model are summarized in Fig. [H 
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Fig. 1. Elementary particles in the standard model. 



The interaction among elementary particles is mathematically described by gauge 
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field theories with the (local) gauge symmetry group of SU{3)c <S) SU{2)l ® f/(l)y, 
where SU{?>)c and SU{2)l ® U(1)y are the gauge groups governing the strong and 
the electro-weak interaction, respectively. The gauge field theory with the SU{?))c 
gauge group is quantum chromodynamics (QCD) [H [2]. Although QCD has been 
established as the underlying theory for the strong interaction, many aspects of it are 
still not well understood, especially in regimes where the strong coupling becomes 
large. The ultimate goal of this work is to improve our understanding about the 
strong interaction in these regimes. Therefore let us begin with a review of properties 
of QCD. 

A. Quantum Chromodynamics 

The fundamental degrees of freedom in QCD are quarks and gluons. Their interaction 
is described by the following Lagrangian, 

Nf 

Cqcd = Y.i'fiil^D, - mf)i^f - \F!i^F;, . (1.1) 

/ 

Here the gluon field strength tensor F^^, reads 

F;, = d.Al - d^Al + igfa,cA\Al , (1.2) 

in terms of the gluon gauge fields A^ (a= 1 ■ ■ • 8). The colored quark fields ^/'j (/ = 
u, d, s, c, b, t) are coupled to the gluons through the gauge covariant derivative 

D, = d,- ig^Al , (1.3) 

where Aq are the Gell-Mann matrices, which are generators of the SU{?))c group 
satisfying 

[KAb\ = fabcK, (1-4) 



fabc being the structure constants of SU{3)c. 

The important properties of QCD include: 

1. Asymptotic freedom. One of the most notable differences between the strong 
and electromagnetic interaction comes from the fact that not only quarks but also 
gluons themselves carry color charge, reflected in the last term of Eq. ( II ■2p . This gives 
rise to gluon self-interactions, which in turn lead to the QCD coupling constant, a^ = 
g'^/4:TT decreasing logarithmically with the momentum transfer in the process, 

.2 ^ 



<Q) = I- 



1.5) 



47r PoHQ'/^Icd) ' 
where Aqcd — 200MeV is introduced as a "non-perturbative" scale where as{Q) for- 
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Fig. 2. Running coupling constant a^ of QCD. Dependence of the QCD coupling con- 
stant, as = g^/An, on the momentum transfer, Q, of the interaction. Figure 
taken from Ref. [3]. 



mally diverges, see Fig. [2J Note that the logarithmic behavior in Eq. f ll.Sp is based on 
perturbative calculations, which break down at momentum transfer well above Aqcd- 
In practice, the scale for the onset of non-perturbative effects is typically given by 



the hadronic mass scale of ~1 GeV. For interactions with large momentum transfer 
Q ^1 GeV (or equivalently at small distance r <^ Ifm according to the uncertainty 
principle), theoretical calculations can be organized in a converging series of terms 
characterized by increasing powers of a^. In this "perturbative" regime, QCD is well 
tested, attaining excellent agreement with experiment. In the opposite limit, due 
to the growing strong coupling constant toward small Q, the perturbative expan- 
sion breaks down and the QCD enters the "strong" regime where "non-perturbative" 
phenomena occur. 

2. Confinement. A striking non-perturbative phenomenon is the confinement 
of color charges, which refers to the fact that isolated colored particles, quarks and 
gluons, have never been observed. They are always confined in "bags" of color-neutral 
baryons or mesons. For example, in deep inelastic scattering (DIS) experiments, when 
a colored quark is knocked out of a proton, qq pairs form along the trajectory of the 
outgoing quark to neutralize color. The knocked-out quark is thus accompanied by 
large concentrations of hadrons along its direction of propagation, so-called "jets" . 

3. Spontaneous Breaking of Chiral Symmetry (SBCS). Usually the u and d (and 
sometimes s) quarks are referred to as "light quarks" ; c, b, t quark are referred to as 
heavy quarks. In the massless limit of light quarks the QCD Lagrangian with two 
light fiavors can be written as 



J^QCD = ^Li^D.^L + ^Ril^D.^n - -FrF^^ , (1.6) 
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J 

4 



where 

/ \ 

u 

^ = ■ (1.7) 

vv 

The iI)l and ipR are the left-handed and right-handed components of the quark field. 
The Lagrangian (II. 6p is unchanged under a rotation in fiavor space {u O d) "inde- 



pendently" with respect to ipi or ipji. This symmetry of the Lagrangian is called the 
chiral symmetry. This symmetry is, however, spontaneously broken by the complex 
structure of QCD vacuum. The latter is filled with various condensates of quark- 
antiquark and gluon fields. In particular, the scalar quark condensate of up and 
down quarks can be quantified by a vacuum expectation value, 

(^^) = {0\tPRtPL + ^lV'/jIO) ^ (-250MeV)3 , (1.8) 

translating into a total pair density of about 4/fm^ (for two flavors). As a result the 
quarks inside the hadrons propagating through the QCD vacuum acquire an effective 
mass, m* ~350MeV, which is much larger than their bare mass, ml^ ^ 5-10 MeV. 
The QCD condensates are thus the main source of the visible mass in the Universe. 
Currently the theoretical efforts of studying QCD in the non-perturbative regime 
are mainly pursued in two directions: 

1. Effective theory/model. Various effective theories are developed with degrees 
of freedom appropriately adapted to specific problems. They are rather successful 
in describing the physics within their applicable energy range. For example in the 
low-energy regime of QCD, chiral perturbation theory was developed with hadrons, 
rather than quarks and gluons, as the effective degrees of freedom. The interaction 
term is dictated by the chiral symmetry, which is approximately respected by the 
original QCD Lagrangian. After the effective coupling constants are determined by 
fitting to experimental data, chiral perturbation theory gains predictive power for 
low energy hadronic reactions/decays. Another example is Non-Relativistic QCD, in 
which the heavy quarks are described by a Schrodinger field theory while the gluons 
and light quarks are modelled by the usual relativistic Lagrangian of QCD. 

2. Lattice QCD. The basic idea of lattice QCD (IQCD) is to study the QCD La- 
grangian in discretized euclidean (imaginary-time) spacetime, with the lattice points. 



called sites, separated by the lattice spacing, a. This effectively introduces an ul- 
traviolet cutoff, A = 1/a on any momentum component. Fermion fields, tp, reside 
on the lattice sites, while the gauge fields. A, are associated with the links joining 
neighboring sites. The IQCD partition function can then be calculated on lattice as 

Z = f [dA][dij][dij]e- ^ ^'^''^ . (1.9) 

Here functional integration denotes summing over all possible field configurations 
with every possible value of the gauge fields, the antifermion fields and the fermion 
fields on each link and sites, respectively. In practice, Monte-Carlo simulations are 
employed to make this summation process possible. In terms of the partition function 
the thermal average values of observables are given by 

{0{A,^,^|J)) = ^j[dA][di,][di,]e-^''''^0{A,^lj,i,) . (1.10) 

Due to the "fermion sign problem" IQCD is currently applicable only in the 
regime of low baryon density. 

As the computer technology and power have advanced in recent years lattice 
QCD has yielded many important results and thus provided profound insights into 
the non-perturbative regime of QCD, see Ref. [4j for a recent review. One of the 
central goals in this work is to establish the link between the IQCD results and 
heavy-ion phenomenology. 

B. QCD Phase Diagram 

A question of fundamental importance is what happens if the hadronic matter is 
compressed so that the distance between hadrons is smaller than the radius of hadrons. 
Intuitively one expects that the boundary of hadrons disappears and the quarks and 
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Fig. 3. Energy density and chiral condensate from lattice QCD. Left panel: Lattice 
calculation of the energy density and pressure of QCD matter as a function of 
temperature, taken from Ref. [5]. Right panel: Strength of the chiral conden- 
sate as a function of temperature, taken from Ref. [6]. 



gluons can move freely inside the entire nuclear matter and become the relevant 
degrees of freedom (deconfinement). It turns out that in the low baryon density region 
this picture is supported by recent IQCD calculations: The energy density shows a 
rapid rise in the temperature region around 170-190 MeV, as shown in the left panel 
of Fig. [3l At high temperatures the energy density e is within 15% of the values 
expected for an ideal gas of quarks and gluons, known as the Stefan-Boltzmann limit, 
which implies the relevant degrees of freedom have indeed transitioned into quarks and 
gluons, forming the quark-gluon plasma (QGP). Lattice QCD calculations also show 
that the rapid increase of e is accompanied by a sudden decrease of quark condensate, 
see the right panel of Fig. [3l implying its evaporation at high temperature similar to 
the evaporation of Cooper pairs above the critical temperature. 

In the regime with low temperature and high baryon density (characterized by 
large baryon chemical potential, /is), theoretical studies [3 E] reveal the existence 
of another deconfined phase, see Fig. HI where the high density quarks form Cooper 
pairs which condense {{qq) 7^0) and result in superconducting of color charge. 
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Fig. 4. Schematic view of QCD phase diagram in terms of the baryon chemical poten- 
tial fiB and temperature T. Figure taken from Ref. [H]. 

Since the formation of deconfined matter (also called quark matter) requires 
either extremely high temperature or density the natural occurrence is expected in 
the early universe or inside compact stars. According to the Big Bang theory the 
universe is believed to have passed through the QGP phase at the age of a few 
microseconds. Inside compact stars such as neutron star the temperature is far lower 
than ISOMeV, but the high density makes it possible that the inside matter is in the 
color-superconducting phase. Since the properties of quark matter are characterized 
by the strong interactions, new insights about the strong interaction is expected to 
be obtained by studying these new states of matter. 

C. Heavy- Ion Collisions 

Since quark matter in nature is far from reach scientists are motivated to create it in 
laboratories, through ultra-relativistic heavy-ion collisions (URHICs). 

The idea is that by accelerating heavy nuclei to very high speed and then colliding 
them a large amount of energy is deposited into a small spatial region and converted 



to thermal energy resulting in extremely high temperature. The currently running 
experimental facilities include the Super-Proton-Synchrotron (SPS) at CERN, Rel- 
ativistic Heavy-Ion Collider (RHIC) at BNL and newly constructed Large Hadron 
Collider (LHC) at CERN. The new Facility for Antiproton and Ion Research (FAIR) 
at GSI (Germany) will be completed within a few years. These experimental facilities 
come with one of two typical setups: the first category is the fixed-target experiments, 
in which only one beam of ions (the projectiles) is accelerated, and its colliding part- 
ner is placed in a stationary target into the path of the beam. The SPS and FAIR 
are fixed target experiments; the second category is the collider, in which both beams 
(projectile and targets) are accelerated and directed to collide with each other. The 
RHIC and LHC fall into this category. Usually the fixed target accelerators have 
higher luminosity (leading to a larger number of collision events per unit time) so 
that more rare reactions can be studied, whereas the advantage of colliders is that 
higher collision energy can be reached (in fixed target accelerators a large amount of 
the energy of the projectile is "wasted" on the kinetic energy of the center of mass of 
two colliding nuclei). 

So far at SPS collisions have been conducted between various ion beams, such as 
proton (p), deuteron (d), O, S, Pb, and different targets such as S, Si, Cu, W, Pb, U, 
at different energies from 20 AGeV to 158 AGeV (for a proton beam it can reach up to 
450 GeV). RHIC has produced collisions between p+p, d-l-Au, Cu+Cu and Au+Au 
at different energies ranging from A/i=22.4AGeV to 200 AGeV. In near future Pb-|-Pb 
collisions will be performed at LHC with up to ■\/i=5.5 ATeV. FAIR will carry out 
heavy-ion collisions with ^/s close to lOAGeV. Different heavy-ion experiments, with 
different beam energies, probe different regions in the QCD phase diagram, Fig.H) the 
matter created in the central region of collisions with higher beam energies is more 
symmetric between baryons and antibaryons, while the lower energy experiments 
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(such as FAIR) enable to study the properties of dense baryonic matter. In this work 
we mainly focus on charmonium production in Pb(158 AGeV)+Pb collisions at SPS 
and ^/s=200 AGeV Au+Au collisions at RHIC. 

The time evolution of a typical heavy-ion collision is sketched in Fig. [51 Two 
Lorentz-contracted nuclei approach each other at close to the speed of light un- 
til primordial nucleon-nucleon collisions occur. After subsequent reinteractions for 
ro=0.5-lfm/c a Quark-Gluon Plasma (QGP) is supposedly created. Driven by the 
pressure gradient the QGP expands and cools (for a duration of tqgp ~3-5fm/c). 
The hadronization then follows with further expansion in the hadronic phase un- 
til the "chemical freeze-out" point when inelastic interactions cease with particle 
abundances fixed; after further expansion/cooling until "kinetic freeze-out" elastic 
interactions stop with particle transverse momentum spectra fixed. The total fireball 
lifetime is approximately 10-15fm/c depending on the beam energy. 
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Fig. 5. Schematic picture of the various stages of a heavy-ion collision. Picture taken 
from Ref. [9]. 

It is convenient to introduce the standard coordinate system for heavy-ion col- 
lisions: The z axis is parallel to the beam line. Since most nucleus-nucleus (A-A) 
collisions are not head-on collisions, there exists a two-dimensional vector connecting 
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Fig. 6. Schematic representation of a noncentral heavy-ion coUision, characterized 
by an almond-shaped initial overlap zone, and a subsequent pressure-driven 
build-up of elhptic flow. Picture taken from Ref. [9]. 

centers of the colliding nuclei in the plane transverse to z axis, which is called the 
impact vector, 6, its length is the impact parameter, h. The x-axis is chosen to be 
parallel to the impact vector, 6, see Fig. [6l The x- and z-axes span the "reaction 
plane" of a given collision. The x- and y-axes span the "transverse plane". The 
component of the 3-momentum of produced particles parallel to 2;-axis is denoted by 
])zi and the transverse component is pf. For relativistic particles it is convenient to 
use the (longitudinal) rapidity instead of the (longitudinal) velocity. The former is 
defined as 

y = tanh"^ (j;) = tanh"^ v, . (1.11) 

Here E = ^Jm? + ff is the energy of a particle. Due to the time-dilation effect parti- 
cles with larger Vz in center of mass frame are "younger" than the ones with smaller 
Vz- This means that the above mentioned evolution of the matter is "measured" by 



the longitudinal proper time T=\/t'^ — z^ rather than by the lab time t, as illustrated 
in Fig. [71 

Due to the short lifetime of the medium, special probes are needed to access the 
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Fig. 7. The spacetime diagram of ultra-relativistic nuclear collisions. Figure taken 
from 1101. 



properties of the medium. The only probes turn out to be the produced particles 
themselves. According to their energies the probes are divided into two categories: 
soft probes and hard probes. 

The soft probes are associated with the particles with relatively low energy (e.g., 
<2GeV), which constitute the bulk medium (>95% at RHIC) created in heavy-ion 
collisions. The soft probes reflect the collective properties of the medium, such as 
thermal and transport properties. One of the key measurements at RHIC is the 
elliptic flow, f2, for the bulk particles {7T,K,p). The elliptic flow, V2, characterizes 
the azimuthal asymmetry of these particles in the transverse plane in terms of the 
second harmonic coefficient of an azimuthal Fourier decomposition of the momentum 
spectra. 



dN 



cPptdy 



dN 



y=0 



TTdpfdy 



[1 + 2v2ipt) cosi2(P) + . . .] . 



:i.l2) 



Here (j) is the azimuthal angle in the transverse plane with (f)=0 for x-axis. At mid- 
rapidity the system is symmetric about y — z plane, so there is no cos term. 

For soft particles {pt <2 GeV) the elliptic flow arises because, in semi-central 
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Fig. 8. Elliptic flow measured by PHENIX and STAR [HI [12] compared to hydrody- 
namic calculations [13] . 



collisions, the geometry of the initial interaction region has the shape of an ellipse, 
see Fig. [61 Once the system thermalizes this initial geometrical anisotropy translates 
into stronger pressure gradients in the direction of the smaller axis of the ellipse. 
This induces momentum correlations among particles which flow preferentially along 
the small axis of the ellipse, leading to a positive V2, see Fig. [8] Since the spatial 
anisotropy is largest at the beginning of the evolution, a measurement of V2 pro- 
vides access to the thermalization time scale tq of the system. Applications of ideal 
relativistic hydrodynamics have shown that the experimentally measured V2{pt) for 
various hadrons (vr, K, p, A) is best described when implementing a thermalization 
time of ro=0.5-l fm/c. 

Other valuable soft probes include particle ratios and the HBT (Hanbury-Brown- 
Twiss) particle interferometry. The former is used to estimate the temperature T and 
chemical potential fiB of "chemical freeze-out" when the inelastic interactions cease 
and the particle ratios are fixed. The latter provides clues about the size, shape and 
time evolution of the medium 
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To further study the properties of the medium, in particular the more microscopic 
aspects, the second category of the probes, namely the hard probes, are required. 

The hard probes are associated with the particles with relatively high energy 
(>2GeV) including either light particles with large momentum or heavy particles 
irrespective of their momenta. Usually a hard probe can only be generated in initial 
hard collisions (their energy scale is usually much larger than the typical temperature 
of the medium) and their initial production can be estimated from p+p collisions. 
By measuring the modification of hard probes after traversing the QGP medium, one 
can obtain information on the microscopic interaction between strongly interacting 
medium and the probe particle. To quantitatively describe the modification due to 
the medium it is convenient to define the nuclear modification factor, 

where A'"coii(&) is number of binary collisions for a given impact parameter, b. The 
relation between N^ow and b will be worked out in Section IVTBIITI The production of 
hard probe (highpj, large mass) particles usually scales with Ncoii{b), implying Raa=^ 
if there are no medium induced modifications. For soft particles Raa is usually less 
than 1 because the production of these particles scales with number of participant 
nucleons, A^part(&)- (^part(&) is smaller than Nco\\{b) in A-A collisions). 

One of the key measurements of hard probes at RHIC is the observation of "jet- 
quenching" |l5j for high pt particles. A jet is a narrow cone of hadrons produced by 
the hadronization of a high momentum parton. If these partons traverse the QGP 
they are expected to undergo coUisional and medium-induced radiational energy loss. 
The energy loss will be reflected in the suppression of high pt hadron multiplicities. 
The data on vr" production in central Au+Au collisions at RHIC have shown a large 
suppression by a factor of 4-5 up to pt — 20 GeV, see Fig. [91 The attenuation of the 
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Fig. 9. RAA{Pt) ratio measured by PHENIX for neutral pions [12], compared to 
jet-quenching calculations [16]. 



parton energy allows one to extract the initial parton density in QGP in terms of the 
transport coefficient q=Q'^/X, which characterizes the (squared) momentum transfer 
per mean free path of the fast parton. 

Another important hard probe is the heavy quark, which is expected to only 
partially thermalize in the medium considering its large mass and the limited fireball 
lifetime. One advantage of heavy quarks over high pt light particles is that they are 
always distinguishable from bulk particles even if they are slowed down. Therefore 
the final heavy-quark (HQ) spectra may be taken as "witnesses" carrying a memory 
of the interaction history throughout the evolving fireball, by operating in between 
the limits of thermalization and free streaming. 

Typical hard probe particles have only one "hard" scale characterized by their 
high energy, therefore they are not particularly sensitive to physics at the energy scale 
of medium temperature, T. However there exists one special hard probe particle which 
has an additional (softer) energy scale (on the order of T) making it very sensitive to 
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physics at the medium temperature. This probe is the charmonium. 

D. Charmonium 

A quarkonium is a bound state of a quark and its own antiquark. A quarkonium made 
of a pair of heavy quarks (c,6) is called heavy quarkonium. Heavy quarkonium in- 
cludes charmonium (cc) and bottomonium (bb). The heaviest toponium does not exist 
because the top quark decays through the electroweak interaction {Tt=l/Tt ~0.1fm/c) 
before a bound state can form. In this dissertation we focus on charmonium which 
can be rather abundantly produced at SPS and RHIC energies. The charmonium 
spectrum in vacuum is summarized in Fig. [TOl Since only vector mesons can couple 
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Fig. 10. Charmonium spectrum. Figure taken from Ref. [T7] . 



to virtual photons and have dilepton as the decay product, which allow for quite 
accurate measurements, in this work we mainly focus on the productions of vec- 
tor charmonium, such as, J/ip and ip'. However we keep in mind that 32% (8%) 
of observed J ftp are from feeddown of Xci'^P') [T8| [19] . which happens at around 
lOOOfm/c, much later than typical thermal medium lifetime (~10fm/c). Throughout 
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this dissertation we will use \E' to denote J/V", Xc and ip' . The vacuum masses for 
J/?/', Xc and Tp' are mj/^=3.1GeV, m;^^Q=3.4GeV, m;^^-^=3.5 GeV, m^^2=3.55GeV 
and m^/=3.7GeV [17]. 

Unlike light quarkonium states the heavy quarks move inside the heavy quarko- 
nium with a speed significantly smaller than the speed of light, with, e.g., (f ^/c^) ~0.25 
for J /ip [20] • As a result the heavy quark bound system can be described with non- 
relativistic Schroedinger approach with static heavy quark potentials, in a similar way 
of describing, e.g., the hydrogen atom. An early ansatz yet successful in describing 
the charmonium spectrum is the Cornell-potential 



4 a 
V(r;T = 0) = --— + ar , (1.14) 

3 r 

with Us —0.35 and a ~1 GeV/fm [23]. The first term corresponds to a Coulombic 
part which originates from one-gluon exchange and is dominant at small distance (r) , 
the second term linear in r reflects the confining interaction. 

As mentioned in the previous section, aside from its large mass, the charmonium 
has another characteristic energy scale, which is its binding energy e^- The charmo- 
nium binding energy is usually counted as the difference between the charmonium 
mass and the open-charm threshold, 

e^ = 2mD — m^ , (1-15) 

with m£, ^1.87 GeV. In vacuum the DD pair is usually considered as the open charm 
threshold for charmonium states. The charmonium states typically have binding 
energies of the order of several hundred MeV, e.g., e^ =640MeV, which is on the the 
same order of typical medium temperatures at URHICs. 

If a charmonium is put inside the deconfined QGP medium the color force be- 
tween c and c is subject to screening by the surrounding colored partons, see Fig. [TTl 
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Fig. 11. Schematic representation of color-Debye screening in a deconfined medium. 

in a way similar to screening of the electric field in dielectric materials: The c{c) 
quark attracts partons in the medium with opposite color charge and forms the "De- 
bye cloud" which screens the color electric field of the c{c) quark. The screening 
effect in the Coulombic part of the Cornell potential can be evaluated with ther- 
mal perturbative QCD (pQCD). In the confining parts it is usually described by a 
phenomenological ansatz in early calculations, leading to the following form of the 
screened Cornell potential at finite temperature [21], 



Voc^ir; T) = -^ (1 - e--(-)0 - ^ e'-^^^ . (1.16) 

A direct consequence of the color Debye-screening is the lowering of charmonium 
binding energies, see Fig. [121 According to thermal pQCD calculations the Debye 
mass is related to the temperature of the medium,T, via 

fin{T) ~ gT. (1.17) 

Here g is the strong coupling constant. Inserting g ~2 (corresponding to a^ ~0.3), 
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Fig. 12. Charnionium binding energies as a function of Debye mass. They are esti- 
mated from Eq. (11.161) . Figure taken from 
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Fig. 13. Free energy of a static cc pair, as computed in lattice QCD [25]. 

we see that above a temperature of T ~ 350MeV the J/ip is not bound any more 
and is expected to dissolve into separate c and c quarks. Based on this mechanism 
J/ip suppression was first suggested in 1986 as a signature of QGP [26j. Similar 
phenomena are expected for the excited charmonium states, such as ip' and Xc- Their 
smaller binding energies imply lower dissociation temperatures. Therefore the entire 
charmonium spectra could provide a "thermometer" for the matter created in heavy- 
ion collisions. 

Recentljo quantitative IQCD computations of the free energy, FQQ{r;T), of a 
static heavy quark pair at finite temperatures, became available [251 [2Z], see Fig. [T31 
This has been used as the main input for recent potential models. However, it remains 
controversial to date whether the free energy or the internal energy, 

dFQQ{r;T) 



UQn{r;T) = FQQ{r;T)-T 



or 



(1.18) 



or any combination thereof, should be identified with a static QQ potential at finite 
temperature, T. Usually the internal energy leads to stronger binding than the free- 
energy. We refer to the former (Vqq = Uqq) and latter (Vqq = Fqq) as strong- and 



^The discussions in this and next paragraph mostly follow Ref. 
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weak-binding scenario, respectively. 

Further progress in thermal IQCD came with the computation of two-point cor- 
relation functions of a quarkonium current, ja, with hadronic quantum number a, 

Ga{r,r) = ((j«(r,f)jt(0,0))) , (1.19) 

as a function of imaginary (euclidean) time, r (also called temporal correlator). The 
imaginary part of the Fourier transform of the correlation function, Ga{T, r), is com- 
monly referred to as the spectral function, 

aa{oj,p) = ImGa{uj,p) , (1.20) 

71 

which is related to the temporal correlator via 

POD 

G^(t,T)= dcoaa{co,T)K{uj,T,T) (1.21) 

Jo 

with the finite-T kernel 

cosh[(c^(r-l/2T)] 
^("^'^'^^= sinh[c./2T] • ^^-22) 

Lattice QCD results for two-point correlation functions are usually normalized to a 
"reconstructed" correlator evaluated with the kernel at temperature T, 

G':%r,T) = / duj a^ico,T*) K{uj,t,T) , (1.23) 

but with a spectral function at low temperature, T*, where no significant medium 
effects are expected. The correlator ratio, 

R^{t, T) = G«(r, T)/G':\t, T) , (1.24) 

is then an indicator of medium effects in Ga{T, T) through deviations from one. Cur- 
rent IQCD calculations find that the correlator ratio, Ra{T,T), in the pseudoscalar 
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Fig. 14. Charmonium correlator ratio in the vector (left panel) and scalar (right panel) 
channel as computed in lattice QCD 



(rjc) and vector {J/ip) channel are close to 1 (within ca. 10%) at temperatures up 
to 2-3 Tc [29| 1301 131] , see the left panel of Fig. [TH In the P-wave channels (scalar 
and axialvector) the correlators ratios are substantially enhanced over 1 at large r, 
see the right panel of Fig. [141 This feature is believed to be due to "zero-mode" 
contributions (at uj=0) which are related to the scattering of a charm (or anti-charm) 
quark, c — > c (or c — > c), rather than to cc bound-state properties [32]. This interpre- 
tation is supported by studies of the r-derivative of P-wave correlator ratios, which 
exhibits a much smaller variation (in the limit that the zero-mode part is a (5-function, 
a^mi^) oc S{ijj), its contribution to the temporal correlator is a constant) [331 [3l]- 

In principle, the in-medium properties of charmonia, such as pole mass, in- 
medium width and dissociation temperature, are fully encoded in their spectral func- 
tion. However, the finite number of data points for the two-point correlator computed 
in IQCD severely hampers the inversion of the transform in Eq. (I1.2ip . rendering the 
determination of the spectral function difficult, although various attempts are made, 
such as the Maximum Entropy Method (MEM) [35] . 

Although associated with large uncertainties, these IQCD data suggest that the 
J/ip bound states might still survive above the critical temperature. As a consequence 
the c and c pairs produced in initial hard collisions may coalesce and regenerate J/ip 
in QGP [361 EZ], rendering the original picture with J/ip suppression as the signal 
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of QGP more complicated. Quantitative calculations disentangling primordial J /ip 
production and regeneration are thus necessary for utilizing J /if) to assess the basic 
properties of the hot and dense medium created in URHICs. These constitute the 
main part of this dissertation. 

E. Outline of Dissertation 

In Chapter [TTl we first give an overview of the physical processes related to charmo- 
nium production in each stage of heavy-ion collisions. Then we introduce our main 
framework: the Boltzmann transport equation approach. The Boltzmann transport 
equation describes the time evolution of the charmonium phase space distribution 
function due to dissociation, regeneration and drifting. If one is only interested in the 
inclusive yield of charmonium the Boltzmann transport equation can be reduced to 
a simpler rate-equation. Next we discuss the important relation between the dissoci- 
ation rate and the regeneration rate, which is the principle of detailed balance. We 
discuss its microscopic origin and show its connection with the equilibrium limit of 
the charmonium phase distribution. We also examine the experimental fact that the 
c and c are always produced in pair and discuss its consequence on the charmonium 
regeneration. Finally we consider the correlation between the c and c in coordinate 
space and estimate its influence on charmonium regeneration. 

In Chapter IIIII we discuss the microscopic interaction between the charmonium 
states and the hot medium. In the QGP we calculate the dissociation rates from a 
"quasifree" mechanism with realistic charmonium in-medium binding energies esti- 
mated from IQCD potentials. We work out the momentum dependence of the dissoci- 
ation rates. Based on the detailed balance we calculate and compare the charmonium 
regeneration rates with different input charm-quark momentum spectra ranging from 
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the hard pQCD spectra to the thermal spectra. With the in-medium width of char- 
monia estimated from "quasifree" dissociation rates we construct a model spectral 
function and calculate the correlator ratios. By comparing these ratios with IQCD 
data we extract the \E' dissociation temperatures, which determine the start time of 
the regeneration process in the kinetic approach. 

In Chapter [IV] we present our description of the space-time evolution of a heavy- 
ion collision, which is modeled by an isentropically and cylindrically expanding fireball 
with input parameters determined from the observed hadro-chemistry and the col- 
lective flow velocities. The resulting thermal evolution scenario is in line with basic 
features of hydrodynamic calculations. In particular we extend the fireball descrip- 
tion to lower energies, which will be explored by future FAIR experiments. In this 
regime it turns out that, despite the low initial temperature, the deconfined phase 
still lasts for about 4fm/c, with a notable reheating process in the mixed phase driven 
by the latent heat. Finally, based on the flow field of the fireball model, we estimate 
the transverse momentum spectra of locally thermalized charmonia. 

In Chapter |V] we estimate the charmonium phase space distribution at the 
moment of thermalization which serves as the initial condition of the Boltzmann 
transport equation. We discuss the primordial production of charmonium in initial 
nucleon-nucleon collisions and its subsequent interaction with the pre-equilibrium 
medium. Emphases are placed on the deviation of primordial charmonium produc- 
tion in nucleus-nucleus (A+A) from scaled p+p due to the so-called cold nuclear 
matter (CNM) effects. 

In Chapter IVII we present the numerical results of our model of charmonium 
production in heavy-ion collisions. We introduce the procedure of applying the ki- 
netic equations to calculate charmonium production in the strong and weak binding 
scenarios. We first compare our results for both inclusive J /ip yields and their trans- 
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verse momentum spectra with data at SPS and RHIC Within current theoretical 
uncertainties we find that both scenarios can reproduce SPS and RHIC data for the 
centrahty dependence of inclusive J /ip production reasonably well. However the par- 
tition of primordial and regeneration yields is quite different in the two scenarios: 
the former dominates for strong binding, while for weak binding regeneration largely 
prevails at RHIC energies (except for peripheral collisions). This difference entails 
the strong binding scenario to be slightly favored by the pt data. We also study the 
effects specifically relevant for high pt J /i' production, including formation time ef- 
fects and i?-meson feeddown. These effects lead to a moderate enhancement of J /ip 
Raa at high pt. Then we explicitly evaluate J /ip regeneration from different input 
charm-quark spectra. Next we compare the results of ip' to if) ratio to experimental 
data at SPS and provide the prediction of Xc and ip' production at RHIC. Last we 
present our prediction for charmonium production at FAIR energies. 

We conclude in Chapter IVIII and give a few directions along which we plan to 
expand the work undertaken in this dissertation. The input charm-quark phase space 
distribution based on Langevin simulations, hydrodynamic simulations of the medium 
evolution as well as a microscopic model for primordial cc and charmonium production 
need to be implemented to improve charmonium toward a more quantitative probe 
of the hot and dense QCD matter. 
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CHAPTER II 

CHARMONIUM TRANSPORT THEORY 
In this chapter we discuss the main framework in which charmonium production in 
heavy-ion colhsions is studied. In Section HTjR] we start with an overview of char- 
monium production in different stages of heavy-ion collisions and relevant produc- 
tion/dissociation mechanisms in each stage. In Section HTJlB] we introduce the Boltz- 
mann transport equation to describe the evolution of charmonium phase space dis- 
tribution functions. The Boltzmann transport equation plays a central role in our 
framework. Finally we investigate constraints imposed by the pair production of 
charm-anticharm quarks and study their influence on charmonium regeneration in 
Section [miCl 

A. Overview of Charmonium Production in Relativistic Heavy-Ion Collisions 

The conventional picture of charmonium {^=J /ip, Xa V"') production in heavy-ion col- 
lisions involves three stages: 1) Hard production stage: The two ultra-relativistically 
moving nuclei collide with each other and charm quark pairs are created from the 
hard collisions between partons. This is a fast process due to the rather large mo- 
mentum transfer (~ 2r7ic=2.5GeV) involved. The typical time scale is estimated 
as r ~ 2^ ~ 0.1 fm/c. 2) Pre-equilibrium stage: If the charm and anti-charm 
quarks produced in initial hard collisions are sufficiently close to each other in phase 
space they can further develop into a charmonium state through non-perturbative 
interactions. The charmonium formation processes are rather slow compared to the 
charm quark production processes. The typical estimates for the charmonium for- 
mation time are ~1 fm/c [20j. In the mean time the partons created in initial hard 
collisions rescatter and are in the process of thermally equilibrating. When the "pre- 
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charmonium" cc pairs travel in this medium some of them may colhde with passing-by 
nucleons and be subsequently dissociated. The process is usually called "nuclear ab- 
sorption" , which will be discussed in more detail in Chapter |Vl The overall duration 
of the pre-equilibrium stage is estimated to be about 0.3 — l.Ofm/c depending on 
the center-of-mass energy of the collision. 3) Equilibrium stage: When the medium 
reaches thermal equilibrium, it is supposed to be either in the deconfined Quark-Gluon 
Plasma (QGP) phase, if the energy deposited in heavy-ion collision is large enough, 
or otherwise in the hadronic gas (HG) phase. The medium continues to expand and 
cool until it becomes so dilute that the particles in it stop interacting with each other. 
This moment is usually called "freeze-out" , which signifies the end of the entire evo- 
lution of the medium. After "freeze-out" the produced particles stream freely to the 
detectors with their yields and transverse momentum spectra fixed at the moment of 
"freeze-out" . The lifetime of the thermal medium depends largely on initial energy 
densities of the medium. Typically it lasts for about 10(12)fm/c for central Pb+Pb 
(Au+Au) collisions at SPS (RHIC). The description of the thermal medium will be 
detailed in Chapter IIVI In the equilibrium stage there exist two competing effects on 
the charmonium abundance: On the one hand charmonia undergo inelastic collisions 
with particles in the thermal bath and are subsequently destroyed, on the other hand 
the charm quarks or D mesons in the medium may recombine and regenerate char- 
monia. Since these two competing processes are inverse processes of each other, their 
relative strength is governed by the principle of detailed balance which ensures that 
the charmonium abundances approach their thermal equilibrium values. The rates 
of approaching equilibrium strongly depend on the temperature and particularly the 
effective degrees of freedom, i.e., the phase, of the medium. In QGP, both the signif- 
icantly reduced ^ binding energies (due to color-Debye screening) and the thermal 
parton density (due to the large color degeneracy) lead to dissociation/regeneration 
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rates much larger than those in HG phase. 

B. Boltzmann Transport Equation 

1. General Setup 

The charmonia in the thermal medium form a system off thermodynamic equilibrium. 
The Boltzmann transport equation is an ideal tool for studying such a system. In this 
work we employ the Boltzmann transport equation to describe the time evolution of 
the charmonia phase space distribution function /^(p, x, t) in the thermal medium [38l 
[391I101I1I1I12], 

dh/dt + v^ ■ V/vt = -a^h + /3* • (2.1) 

Here, the distribution function fq, (x, p, t) is the number of charmonia (\1/ = J /'ip^ Xc, V"') 
per unit phase space volume at a given time t, 

^*^"'^'^) = ^%d^- (2.2) 

The first term on the left-hand side of Eq. (12. ip describes the rate of change of 
the charmonium distribution function at a specific position in the phase-space. The 
second (drift) term reflects the charmonium diffusion from one position to another 



with velocity v^^ = p / yjp^+m^ . Since the mass of \l/ is much larger than typical 
medium temperature (m^ ^ T) we neglect the elastic scattering between \1/ and the 
medium and assume the \1/ to move on a straight line. Therefore the left-hand side 
of Eq. (12. ip represents the change of charmonium density other than charmonium 
diffusion. Other than diffusion the change in \l/ distribution function can only be due 
to the dissociation and regeneration of charmonia, which are accounted for by the first 
( "loss" ) and second ( "gain" ) term on the right-hand side, respectively. The loss term 
consists of the product of the charmonium phase space density and its dissociation 
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rate, aq,{p,x). The gain term, Pq,[p,x), reflects the charmonium regeneration rate 
from charm quark coalescence. 

One of the premises of Boltzmann transport equation is that the medium should 
be "dilute" enough so that a subsequent collision occurs well after the end of the first 
one (no quantum interference between successive collisions). The typical duration for 



one collision can be estimated by the reciprocal of the virtuality, a/Q^ — m'^, of the 
^ after scattering with thermal particles, which can be estimated as follows: Denote 
the four momentum of the \E' before a collision as (p^/(2r?7„j,) + m,j,,p), and the four 
momentum of a thermal parton as {k, k) with k=\k\. Then the four momentum of ^ 
after scattering becomes (go, g) = (p^/(2m,j) + m^ + k,p + k). The virtuality is given 

by 

Q'^ -m% = ql- ^ -m% = ( V rrm, + k] - {p + kf - m% 



p^k p"^ 



2m^k -2p-k + ^ + -^ . (2.3) 



Assuming thermal momenta for \E' and the parton, so that p^ /{2m) r^ T, p r^ y/m^T, 



k ^ T, the virtuality y/Q^ — m\ to leading order in m,^(^ T) is on the order of 
\/m<i,T. Therefore with a typical charmonium mass ~3 GeV and a typical tempera- 
ture T ~0.3GeV, the estimated duration for a dissociation process is around 0.2fm/c. 
On the other hand the typical \E' dissociation rate at RHIC energy is on the order 
of lOOMeV, as will be evaluated in Section IIII||B| corresponding to an average time 
of 2fm/c between successive collisions. Therefore the average interval between \E' in- 
elastic collisions is indeed much longer than typical collision time thus satisfying the 
premise of Boltzmann transport equation. 

We proceed to work out the microscopic expression for the dissociation and gain 
rates, a<^ and /3^. From now on we will omit the vector sign "" " on top of x 
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and p if no ambiguity arises. For simplicity let us first consider a 2-body dissocia- 
tion/regeneration process, e.g., \E' + (7 ^ c + c, in which a \E' scatters with a gluon 
and is subsequently dissociated into a cc pair. The dissociation rate is expressed as 
follows: 

a<i,{p,x) = ^^ / d^i{pg) d^2{Pc,Pc} 

X (27r)^5(")(Pc +P-C-P- Pg) dg fgiPg) \M^g^ccis, t) | ^ , (2.4) 

where E\^ = ^Jp^ + ?7i|, is the charmonium energy, 

is the gluon phase space distribution function with dg=\^ being the color-spin degen- 
eracy factor for gluons, and lA^cc-s-^gl^, t)p is the transition matrix element between 
the initial state, \E' + (/, and the final state, c ^ c, as a function of s = (p + PgY 
and t = {p — Pc)"^. The "bar" on top of |A^p stands for summing over color and 
spin degeneracies of final states and averaging over color and spin degeneracy of 



initial states, e.g., |A^^g_s.ccP=S|A^*g^ccP/('^g'^*), with dg and (i^ being color and 
spin degeneracies of the gluon and the charmonium. The UlAix^jg^csl"^ denotes the 
transition matrix element with both initial and final state degeneracy summed over. 
In Eq. (12.41) d^i{pg) and d^2{Pc,Pc) are Lorentz invariant 1-body and 2-body phase 
space integration measures, 

d T) d J)- 

^^^(^-^^")=(2.)32E.(2vr)32E,' ^'"^^ 

where Eg, Ec and Ec are the gluon, c and c energies, respectively. 

It is convenient to express the dissociation rate a.^ in terms of the charmonium 
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dissociation cross section a^, which can be expressed in terms of the transition prob- 
abihty, |A^ cc^^gp, as follows, 

If 1 

a^,{s) = —— / rf$2(pc,Pc)(27r)V^)(p,+p5-p-p<,)77^ :\M^g^c-c{s,t)\'^ , 

/-C/Xl; J Zl^g\V\^ — Vg\ 

(2.8) 



where Vx^, = p/\/j^+rn^ and Vg = Pg/ y^p^'+rn?^ axe the velocities of the charmonium 
and gluon, respectively. By comparing Eq. fl2.8p with Eq. fl2.4p one can identify the 
relation between a^ and (T<^, 

/d T) 
7^^dgfg{pg)Vrel(T'S,is), (2.9) 

where Vrei=\viii — Vg\ is the relative velocity between the gluon and the charmonium. 
Next we proceed to the microscopic expression of the gain rate, P^, given by 

/3*(p,a;) = ^^ / d^i{pg)d^2{Pc,Pc){'2n)^S^^\pc + Pc - P - 

X \M,,^^g{s,t)\^ 4 Mp,) 4 /c-(Pc-) (1 + fgiPg)) , (2.10) 

where dcfdPc) and dcfc{Pc) are charm and anti-charm phase space distribution func- 
tions, 

with dc=dc=Q being the color-spin degeneracy factors for c and c, respectively. The 
factor of 1 + fg{Pg) is the Bose-enhancement factor for the final state gluon. The 



"bar" on top of |A^cc^*g| denotes summing over final states and averaging over the 



initial states, \M.cc^^ig\=Yl\M.cc^<i,g\/{dcdc). We note that the underlying dynamics 
of the transition between ^'-l-^E' and c + c^ QCD, is time-reversal symmetric, therefore 
the transition probability satisfies the principle of detailed balance as 

S|A^,e^^,|2 _ Y.\M^g^cc? . (2.12) 
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It is straightforward to extend the above expressions for the dissociation and gain 
rates to 2-to-3 processes, z + \I/— j-z + c + c {i=g,q,q) and its inverse 3-to-2 processes, 
i + c+c ^ i + \l/ . What needs to be modified are: 1) Replacing the 2-body integration 
measures in Eqs. (]2.4p . (J2.8P and (I2.10p by 3-body integration measures, 

d^Pc d?Pc d?Vi 

dUp.,P,,P.) = (2^)32,,^ (2vr)32E, (2vr)32E. ' ^^''^^ 

accounting for the phase space of the extra hght parton in the final state of dissociation 
(or the initial state of regeneration), 2) Inserting the four-momentum of the extra light 
parton into the 4-momentum conserving delta function. 3) Supplying the phase space 
distribution function, difi{pi), for the extra parton if it appears in the initial state 
or the Bose-enhancement/Pauli-blocking factors, 1 ± fi{pi), if it appears in the final 
state. The "±" takes +(— ) ior i = g {i = q,q). The exphcit expressions of a^, a,^ 
and /3,j, for two-to-three processes are 

(^^(s) = J^Yl / d^3{Pc,Pc,Pi)i27r)^5^^\pc + Pc + Pt - P - Pi) 



X — -\M^,^Ms,t)ni±Mp^)) , (2.14) 

Zll/g\V\^ Vi I 

a^{p,x) = ^^^ / d^i{pi)^3{pc,Pc,pi){27i)'^6^^\pc+Pc + Pi-p-Pi) 



X \M^i_,-4s,t)\^dJ,{p,){l±f,{Pi)) , (2.15) 

f3^{p,x) = ^^J2j c?$i(p»)c/$3(Pc,Pc-,Pi)(27r)V^)(p, + p5 + . 



Pi-p-Pi) 



X \Mcs^<i^i{s,t)\^ dc fc{Pc) dc fc{Pc) di fi{pi) {1 ± fi{pi)) . (2.16) 

Here pi denotes the momentum of the light parton for the three-body state (together 
with c and c). 

We will illustrate in Section [ITTIIC] that in the case where the charmonium binding 
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energies are small compared to the medium temperature, a 3-to-2 regeneration pro- 
cess, -i + c+c — )■ z-|-\l/ (z = g,q,q), can be readily factorized into a 2-to-2 (perturbative) 
quasi-elastic scattering process, i + c{c) — )■ i + c'{c') and 2-to-l (non-perturbative) co- 
alescence process, c+c — ?■ \E' owing to the different time scales for these two processes. 
The four-momentum is conserved for the entire 3-to-2 regeneration process. 

So far we have established the connection between the charmonium transport 
equation and the microscopic dynamics of \1' dissociation/regeneration. The required 
inputs from microscopic calculations are 1) The transition matrix element, | A^jcc-i-'i'iP; 
2) Charmonium binding energies, e$, which determine the initial state phase space 
for the dissociation processes and the final state phase space for the regeneration 
processes; 3) The in-medium charm quark spectra, fdPc), fc{Pc), determining the 
charmonium regeneration rates. We will discuss these input quantities based on mi- 
croscopic calculations in Chapter UTTI With all these quantities supplied by microscopic 
calculations we are ready to solve the Boltzmann transport equation (12. ip . 

The Boltzmann transport equation, Eq. (12. ip . is a first-order, linear partial differ- 
ential equation. It can be solved using the method of change of variables. Introducing 
the new variable 

u = X — v^it ; (2.17) 

we have 

U{x,t) = U{u + v^t,t) , (2.18) 

a*(a;,t) = a^(u -I- f^t, t) , (2.19) 

I3m{x, t) = I5m{u + vmt, t) . (2.20) 

Since in Eq. (12. ip there is no derivative with respect to p, we consider p as a parameter 
and suppress its dependence in /^, a^, /3^ for the moment. Substituting Eq. (I2.17P 
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into Eq. f l2.ip . we obtain 

^f^iJ{u + vq,t,t)/dt = —aii,{u + Vii,t,t)fq,{u + vqit,t) + f3q,{u + vqjt,t) . (2.21) 

Now Eq. (I2.2ip is a first-order ordinary differential equation with u as a parameter. 
It has the following solution, 

Mu + v^t,t) = /*(w + t;^to,to)e-^'^'*'"*^"+^**''*'^ 
i-t 
+ dt%(M + i;*t',t')e"^''''*""*^"+'**"'*"^ (2.22) 

Jto 

Finally we substitute Eq. f l2.17p into Eq. f l2.22p and obtain the solution for the Boltz- 
mann transport equation. 

Mp, X, t) = Mp, x-v^{t- to), to)e-^o ^*'-*fe-^*(*-*')/) 

+ / dt'l3^{p,x- v^{t - t'),t')e-/.'^*""*(p.--*(*-*")/') , (2.23) 

where to is the start time of the evolution. The initial phase space distribution 
of charmonia, /^(p, x,to); is determined by initial hard production and cold nuclear 
matter effects, which are the main topic of Chapter |Vl The solution of the Boltzmann 
transport equation, Eq. ( I2.23p . consists of two terms: The first term decays with time 
and describes the dissociation process of the initially produced charmonia. This 
component of the solution is usually referred to as the "primordial" component. The 
second term increases with time and describes the regeneration process of charmonia 
from coalescence of charm quarks in the medium. This component is referred to as 
the "regeneration" component. 
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2. Rate-Equation 

If one is only interested in the inclusive yield of charmonia, Ni^,, it is convenient to 
employ a rate-equation instead of the more differential Boltzmann equation. In this 
section we derive the widely employed rate-equation from the Boltzmann equation. 

We start by integrating over the entire charmonium phase space on both sides 
of the Boltzmann equation fl2.ip . The left-hand side becomes 



727r)3 '^^•^**^^' ^' ^''/^^ + ^* ■ '^•^*'^^' ^' ^)) 

^ f j3^ d^P -c /„ ^ -t\ , f d^P 



j d'x-0^Mp,x,t) + j -0^ j d^xV ■{Mp,x,t)v^) 



= dN^,{t)/dt+ / — -^ (h dS ■ {fq,{p,x,t)v^) 

= dN^,{t)/dt . (2.24) 

In the second equality we used the definition of charmonium phase space distribution 
function, Eq. (12. 2 p and Gauss's law, where S" is a large surface enclosing the inte- 
gration volume of coordinate space. In the third equality we used the fact that the 
charmonium phase space distribution fq,{p,x,t) drops sufficiently fast as |x| — )■ oo. 
Next we proceed to the right-hand side, which consists of the loss and the gain terms. 
For the loss term, if spatial homogeneity is assumed, namely, a^{p,x,t) — )• a^(p, t), 
we have 

d V f d V f 

d^x——^aq,{p,x,t)fq,{p,x,t)= / ——a^{p,t) / d^x fq,{p,X,t) 

r d T) 
Ria^{{p),t) / -—— d^xU{p,x,t) 



J (27r)3 
^r^((p),t)iV^(t), (2.25) 

where (p) is the average momentum of \E', and the aii,{{p),t) is conventionally denoted 
as r^((p),t). Similar manipulations can be apphed to the gain term. With the 
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assumption of spatial homogeneity the integration over coordinate space reduces to 
the multiphcation with the volume of the medium (fireball), Vfb- The integration 
over charmonium momentum can also be performed. The resulting integrated gain 
term is conventionally denoted by G<ij{t), 

r d T) 

G'^(t) = / d^x (3q,{p,x,t) 



(2vr) 



= VpBJ^,P^{p.t). (2.26) 

Putting Eqs. ([221!), '12^241) and (12:261) together we obtain the rate equation [13] de- 
scribing the time-evolution of the inclusive charmonium yield, 

^^^ = -T^{t)N^{t) + G^{t) . (2.27) 

dt 

3. Equilibrium Limit 

Let us consider the equilibrium limit of the solution of the Boltzmann transport 
equation, Eq. (12. ip . which is defined by 

mp) = ^ , (2.28) 

a^{p) 

where f^{p), P^aij)) and a^(p) are also assumed to be homogeneous in space. It is 
easy to verify that the equilibrium distribution f^{p) solves the Boltzmann equation. 
In this limit the (integrated) gain term in the rate-equation, Gq,, can be written in 
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terms of f^{p) as 

/dp 

= T^{t)Nf{t) . (2.29) 

Substituting Eq. fl2.29p into Eq. f l2.27p . we obtain another common form of the rate- 
equation [S], 

^^^ = -T^{t)[N^{t) - N'^{t)] . (2.30) 

This form is particularly convenient in the special case where the momentum spectra 
for both the charm-quark and the light partons are thermal. In this case the inclusive 
yield of charmonium in the equilibrium limit, f^{p), is independent of the transition 
probability, |A^p. It is instructive to further discuss this special case. 

For simplicity we first consider a 2-to-2 process, "^ + g ^ c + c. We denote 
thermal (Boltzmann) charm-quark spectra as 



fciPc) = 7c e- V^^^^/^ = ^, e-^^/^ , (2.31) 

where 7c is a charm quark fugacity reflecting its total yield deviating from chemi- 
cal equilibrium (since rric ^ T we cannot expect charm quarks to reach chemical 
equilibrium). The gluons follow the thermal Bose distribution. 



exp ^^f^ - 1 exp I - 
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where Eg is the energy of the gluon. Plugging these spectra into Eq. fl2.10p . we obtain 

P^{P) = ^^ / d<^i{pg)d<l>2{Pc,Pc){27c)^S\pc + Pc - P - . 



X \Mcc^^g)\^dcfc{Pc)dcfc{Pc) (1 + fgiPg)) 
ig I t/$i(p,)rf$2(Pe,Pc-)(2vr) V(pe + P, ' P ' Pg) 

:^ I d^,{pg)d^2{Pc,P,){27ry6\p, + p,-p-pg) 
y< \M—^g\' e-^^--^^^^^ (1 + fgiPg)) 



2Eq, 



(i$l(Pg)d$2(Pc,Pc-)(27r)4r(Pe + Pc-P-Pg) 



y<\M^g^,,\'e~V<^/^dgfg{Pg) 

= -f^,d^e-V<+py^ a^{p) , (2.33) 

where in the third and fourth equahty we used four-momentum conservation im- 
posed by the 4-D (5-function, and an identity for the Bose-distribution for the gluon, 
g-^s/^^l _l_ fg{j)g)) = fg{pg). In the fifth equality the detailed balance between 
Ti\M.^,g^cc\'^ and S|A^cc-i.'^gP, Eq- (I2.12p . and the definition of dissociation rate a^, 
Eq. (12 ■4p . are used, The color-spin degeneracy factors for \1', gluon, c and c are (i^=3, 
(ig=16, dc=dc=Q. Although this derivation is made for 2-to-2 processes, it is straight- 
forward to verify that this relation also holds for 2-to-3 processes as long as the 
detailed balance and four-momentum conservation hold. Comparing Eq. (12.331) with 
Eq. (I22HD, we obtain 



mp) = l!d^e-V<^/^ . (2.34) 

This relation shows that as long as both charm quarks and light partons are in 
thermal equilibrium the regenerated charmonia have the thermal momentum spectra 



39 



as their equilibrium limit which does not depend on the microscopic details of the 
dissociation/regeneration mechanisms. Different microscopic transition mechanisms 
will only affect the transition rates, namely the speed with which charmonium spectra 
approach the equilibrium limit. 

Since the thermal production and annihilation rates of cc are believed to be 
small in heavy-ion collisions at SPS and RHIC energies, cc pairs are assumed to 
be exclusively produced in primordial N+N collisions and conserved thereafter [36] . 
The charmonium equilibrium limit at a given time can be conveniently evaluated by 
solving the following charm conservation relation for 7c, 



V 2 J (27r)3^'"'^"^ '"'^"^^ J (27r)3- 

= 7cnf + i!nf , (2.35) 

where Ncs is the total number of cc pairs from primordial production and, n^^ and 
n^ are thermal spatial densities for c and \E', 



th J / d.^Pc V"'c+P% 



< = dc I ^e ^ , (2.36) 



th , f d^P V'"*+''" 

* " * / W)"^ — ^~ • ^^■^'^^ 

This relation shows that in thermal equilibrium the partition of all charm quarks be- 
tween the open-charm and charmonium states is solely determined by their respective 
masses and degeneracy factors |15] . 

We should keep in mind that the charmonium equilibrium limits have this nice 
feature only if the charm quark spectra are thermal, which, however, may not be the 
case in heavy-ion collisions, since charm quarks are heavy and take a rather long time 
to thermalize, compared to the typical lifetime of the medium (delayed by a factor of 
rric/T compared to light partons). We will discuss more realistic charm quark spectra 
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and their implications for charmonium production in the medium in Chapter IIIII 

C c-c Correlations for Charmonium Production 

1. Charmonium Production in the Canonical Ensemble 

Before we finish this chapter and discuss the input to the kinetic equations from micro- 
scopic calculations we need to address another issue associated with '^ regeneration. 
The problem originates from a simple experimental fact: The charm and anti-charm 
quarks are always produced in pairs. Its impact on charmonium production can be 
easily seen in a situation where few charm quark pairs are produced per event. For 
example, at SPS energy there is one charm quark pair produced every ten events 
in central Pb-|-Pb collisions, namely on average Nc=Nc=Q.l per event. Accordingly 
one expects the average number of regenerated ^ to be 1/10 of what we would get 
if there is one charm quark pair created in each event. However Eq. fl2.10p naively 
suggests that the average number of regenerated \E' is proportional to NcNc=Q.Ql, one 
hundredth of what we would get for the Nc=Nc=l case, an underestimate by a factor 
of 10! This example illustrates the importance of the correlation between charm and 
ant i- charm production in charmonium regeneration. 

To systematically solve this issue we employ the statistical description of charm 
quark pair production in a grand-canonical ensemble, where all thermodynamic prop- 
erties of the charm pair system can be derived from its grand-partition function Z |46j . 
Let us start with the 1-body partition function (for 1 charm or anti-charm quark), 

^1 = IcuTVfb = Ic Vfb 4 / ^e- V^^^^/^ , (2.38) 



J (2vr)3 

where VpB and T are the volume and temperature of the system, the charm quark 
fugacity 7c is to account for the chemical off-equilibrium of charm quarks. K2[x) is 
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the modified Bessel function, n^^ is tlie tfiermal density of charm quarks, defined in 
Eq. f l2.36p . The partition function with k charm quarks is 

^fe = If , (2.39) 

where k\ in the denominator resuhs from the indistinguishabihty of k charm quarks. 
If we now impose the constraint that for each charm quark there exists its partner - 
anti-charm quark, we obtain the partition function for k charm pairs, 

Zr = Z, ■Zk = ^^. (2.40) 

The k (anti-) charm quarks are indistinguishable from each other, however the charm 
quarks are distinct from anti-charm quarks. 

Next we sum over the contributions from arbitrary k pairs of charm quarks and 
obtain the grand-partition function Z, satisfying the constraint N^ — Nc=0 for each 
individual event. 

Z = J2zr = j:^ = WZ,), (2.41) 

fc=0 fc=0 

where Iq is the modified Bessel function. Sometimes such an ensemble is called 
"canonical ensemble" . We keep in mind however that this canonicality is only referred 
to the strict constraint of the net charm number (A^^ ~ ^c =0) rather than the "total" 
charm number, Nc + N^, which can still fluctuate event by event. 

With grand-partition function obtained we are ready to evaluate average values 
over events (denoted by "(■■■)" in this section) of any thermodynamic quantity as 

OO 

(0)^^0(A:)P(A:), (2.42) 

fc=0 

where 

ypair 

P{k) = ^— , (2.43) 
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Fig. 15. The probability distribution of integer number (fc) pairs of charm quarks 
in one event. Left paneh {Ncc)=0-l corresponding to central Pb+Pb col- 
lisions (within Ay=1.8 around y=0) at SPS {y/s=17 .3AGeV) . Right panel: 
{Ncc)=5.0, corresponding to central Au+Au collisions (within Ay=1.8 around 
y=0) at RHIC (v^=200AGeV). 

is the probability of an event with an integer (k) number of pairs of charm quarks. 
For example, the average number of (open) charm pairs is 

k=0 

Matching (Nop) to the experimentally measured (Ncc) allows the determination of Zi, 
Eq. (I2.38P and therefore P{k) (the number of the charmonium states in the system 
is numerically negligible compared to (Nop)). 

In central Pb+Pb collisions at SPS {^/s=17 .3AGeV) and in central Au+Au col- 
lisions at RHIC {^/s=200AGeV) there are on average 0.1 and 5.0 charm-quark pairs 
produced in the rapidity window of Ay =1.8 around y=0 produced, respectively. The 
corresponding probability distributions P{k) are shown in Fig. [151 

Eq. (I2.23P states that the number of regenerated \E' is proportional to the gain 
rate, /9^, which, in turn, is proportional to the product of the charm and anti-charm 
quark phase space distribution, fcfs- In events with exactly k pairs of charm quarks 
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produced, the regeneration component can be written as N^^{k) = k'^Nl'^^, with A^^*^^ 
being the number of \l/ regenerated from 1 cc pair. Therefore the average number of 
charmonia regenerated in the canonical ensemble over k is 

{ND = Yl N'^''ik)Pik) = Yj^Yl ''"^r = K'Z! . (2.45) 

Specifically, in the equilibrium limit we have 

Nr^ = NP = 7ji)V>,. / (1^ d^e-V^<^/^ = 7J1) VpB^^ , (2.46) 

where 7c(i) satisfies 

l-7.„^.B/|^4e-V^'-^BlI^. (2.47) 

Substituting Eq. (12:471) into Eqs. fl2:i6|) and (jMS]) we obtain 

(iV;"^) = N^Z', = ^IVfbu^ . (2.48) 

With {Nop) and {N'^) known we are ready to write down the charm conservation 
equation in the "canonical" ensemble [17], 

= IcncVpE-Tf 7^^ + Ic^FBns^ ■ (2.49) 

This equation allows us to solve for 7c and obtain (N^). Again the second term on 
the r.h.s is numerically negligible due to m^ ^ rric- 

It is instructive to examine two limits where a large (small) number of charm- 
quark pairs are produced. We first note the following property of Bessel functions. 



(2.50) 
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Therefore, in these two hmits, Eq. f l2.44p reduces to 



Zx , {N,^) > 1 
{N,,)^{ (2.51) 

Z^ , {N,,} < 1 . 

Plugging Eq. f l2.5ip into Eq. (12.45 p . one obtains 

[ {N,,)^Nr^ , (N,,) > 1 
{Nr)^{ (2.52) 

In the hmit of (Ncs) ^1, the yield of the regeneration component is proportional 
to (Ncc)"^, the same result one would get if the canonical constraint Nc = Nc was 
neglected; therefore the this limit is often called the "grand-canonical" limit. In the 
opposite ("canonical") limit, where (Ncs) <^1, the regeneration component is propor- 
tional to (Ncc) '■ the canonical constraint Nc = N^ effectively enhances the charmonium 
regeneration by a factor of l/(A^cc) over the grand-canonical limit. As a side remark, 
particle production in the "canonical ensemble" can also be conveniently described 
with a kinetic master equation approach as developed in Ref. [IBj, which has been 
applied to study strange particle production in low energy heavy-ion collisions |1U] 
where the net "strangeness", {Ng — Ns), is constrained to zero in each event. 

2. Charm-Quark Correlation Volume 

It turns out that the correlation between the charm and anti-charm quarks goes 
beyond "production- in-pair" in each event. The charm and anti-charm quarks are 
produced at the same spatial point in hard N-l-N collisions and then recede from 
each other. Before the medium evolution stops (freeze-out) they can only diffuse 
into a limited portion of the fireball volume. The volume they explore is referred 
to as charm-quark correlation volume [H], Vco, which in general is smaller than the 



45 



full fireball volume, VpB- A schematic estimate of Vco(t) will be presented in Sec- 
tion IVIlRI The restriction of cc pairs to within the correlation volume effectively 
increases the cc coalescence probability: One charm quark can more easily find its 
partner nearby inside the correlation volume Vco with the effective spatial density, 



n 



eff _ 1 f d^p reff 



(^c! T^^ft {v)i of its partner given by 



c — "-c J (2tt) 



"=" -w,ri' P-^^> 



'CO "CO 



for an event in which k charm quark pairs are produced. The effective density, nl , 
is larger than k/VpE (the c density without correlation volume effect) by a factor 
of VpB/{kVco)-, if we assume k anti-charm quarks residing in k correlation volume 
"bubbles" (no merging of correlation volumes). 

It is instructive to work out the parametric dependence of the regeneration com- 
ponent, N'^^^lk), on the charm quark pair number k in two limits: 

(a) In the limit of {Ncc) ^ Vfb/Vco the overlap of correlation "bubbles" can be 
neglected; we have N'^f^ik) = [VpE/ikV^o)] ■ k^Nl''^ . 

(b) In the opposite limit, (A^cc) ^ ypB/yco, the k correlation volume "bubbles" 
maximally overlap with each other with their total volume filling the entire fireball 
volume VpB, leading to an effective anti-charm quark density n^^ = k/Vps- In this 
limit there is no correlation volume effect so we have N^'^^{k) = k'^N^^^. In these two 
limits the average number of regenerated ^ is 

Comparing Eq. f l2.54p with Eq. fl2.44p one sees that in the case of 1 ^ {Ncc) ^ 
Vfb/Vco the correlation volume effect renders (N^^) to depend linearly on (Ncc) even 
if the cc system is in the "grand-canonical" limit. In other words, the correlation 
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volume effect "delays" the transition between the "canonical" and "grand-canonical" 
ensembles from (A^ce) ~ 1 to (A^cc) ~ Vfb/Vco- In the general case where (A^cc) is 
comparable to Vfb/Ko the merge of correlation volumes needs to be taken care of on a 
term-by-term basis in the series of N^^^{k). Whenever two correlation volume bubbles 
merge their new total volume (smaller than the sum of their individual volumes) 
should be used for estimating the effective anti-charm density n^-'^'^. In Section IVTIOTI 
we will study the sensitivity of J/tp regeneration on different correlation volume sizes 
based on a rather schematic prescription of merging correlation volumes: If at any 
given time, t, the total correlation volume of individual "bubbles" kVco{t) is larger 
than the fireball volume VpBit), kVco{t) is set to Vpsit) for determination of the 
effective anti-charm density n^-'^'^ in the subsequent evolution. 

We conclude this section by pointing out that a more rigorous and systematic 
way to account for canonical-ensemble and correlation volume effects is to use the 
joint phase space distribution function, (/c/g), as the input for \E' regeneration rate, 
/3^. The (fcfc) can be obtained, e.g., from high statistics Langevin simulations of the 
evolution of charm quark pairs in the medium. 
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CHAPTER III 

CHARMONIUM IN THE HOT MEDIUM 
In this chapter we discuss the microscopic mechanisms for the reactions of charmo- 
nia, \& {^=J /ip, Xci '4'')i i^ the hot medium created in heavy-ion collisions. In Sec- 
tion lHIlAl we review the relevant properties of \I' in the QGP relevant for charmonium 
production. In Section [Inl[B| we discuss the charmonium dissociation mechanisms for 
^! with reduced in-medium binding energy tB < T, the quasifree-dissociation pro- 
cess of charmonium dissociation, and compare to the more traditional mechanism, 
gluo-dissociation process. In Section UTTimi we explicitly calculate charmonium regen- 
eration rates in QGP with several input charm-quark spectra. In Section IIII|[D| we 
briefly discuss charmonium dissociation in the hadronic matter (at T below the crit- 
ical temperature Tc). In Section ITnllEl we construct charmonium spectral functions 
based on the quasifree dissociation rates and compare to lattice calculations, through 
which we extract charmonium dissociation temperatures, T^^^. 

A. Charmonium in QGP 

On the microscopic level the medium affects charmonium states in three ways: 

1) Debye screening: The QGP is a deconfined medium therefore the freely moving 
colored partons screen the binding force of a charmonium state and decrease its 
binding energy, esiT). The in-medium esiT) is defined as 

eB(T) = 2ml{T)-m^(T) , (3.1) 

where m* and rnqt are in-medium charm-quark and charmonium masses. The in- 
medium part of m* is identified with the asymptotic value of the heavy-quark (HQ) 
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Fig. 16. Temperature dependence of in-medium charm quark mass in the strong (sohd 
hne) and weak binding (dashed hue) scenarios. Figure taken from Ref. [28] . 



potential, as displayed in Fig. [161 



m:{T)^ml + VQQ{r^oo-T)/2. 



(3.2) 



In this work we approximate the HQ potential, VQQ{r\ T), based on the IQCD heavy- 
quark free energy, FQQ{r;T), (recall Section UllDll within either the "strong" {Vqq = 
Ugqir^T)) or "weak" binding {Vqq = FQQ{r;T)) scenario. The in-medium masses 
decrease with temperature appreciably, while the magnitude of m*{T) is significantly 
smaller in the weak-binding compared to the strong-binding scenario. 

For a quantitative estimate of esiT) we take recourse to the potential model of 
Ref. [50] where quarkonium spectral functions and correlators (recall Eq. (11.201) and 
Eq. ( Il.igp ) have been calculated in a thermodynamic T- matrix approach, consistent 
with vacuum spectroscopy and including relativistic corrections for a proper descrip- 
tion of scattering states. The calculations in there have been carried out for both 
free and internal energies as potential, and for two different IQCD inputs [271 El]- 
Since the internal energy leads to stronger binding than the free energy we refer to 
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Fig. 17. Temperature dependence of J /ip binding energy in the strong (solid line) and 
weak binding (dashed line) scenarios. 

the former and latter as strong- and weak-binding scenario, respectively. In both 
scenarios (and for both potentials), an approximate constancy (within ±15%) of the 
correlator ratios, Eq. (II. 24^ . for pseudoscalar charmonium has been found (see lower 
panels of Fig. 12 and 14 in Ref. [50]). We believe that these results provide a rea- 
sonable bracket for potential- model results. Similar to m*{T), the resulting binding 
energies (plotted in Fig. [T7j) also decrease with T, again being significantly smaller 
in the weak-binding scenario. These features are, in fact, the main reason that both 
scenarios can be compatible with the small variations found in the IQCD correla- 
tor ratios: for weak/strong binding, a small/large constituent mass combines with a 
small/large binding energy, respectively, leading to an approximate compensation in 
the bound-state mass, niiiiiT), recall Eq. ( 13. ip . 

2) Collisional dissociation/regeneration: Charmonia are subject to dissociation 
through collisions with particles in the medium even if they have finite binding en- 
ergies. We assume both chemical and thermal equilibrium for light partons, so that 
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their pt spectra follow the thermal Bose/Fermi distribution, 

fhik,T) = A_ , t = g, (3.3) 

exp (^^ j - 1 

fp{k,T) = /^2-^^ ' i = u,u,d,d,s,s . (3.4) 

exp ' — — 



Their thermal masses rrii are guided by perturbative QCD (pQCD) calculations, 

^l = ^1 + -^r- ' (3-6) 



6 



r-l = (1 + ^)^, (3-7) 



where Nf=2.5 is the number of flavors, and the coupling constant g is adjusted so 
that the resulting energy density of the QGP medium, e(T), reproduces the lattice- 
QCD results, see Chapter |IV] for more details. The microscopic mechanisms of the 
scattering between "$ and particles in the heatbath will be detailed in the next section. 
3) Bose-enhancement/Pauli blocking: In charmonium dissociation or regener- 
ation processes light partons may be produced in the final state, e.g., i + \E' — )■ 
i + c + c{i = g,q,q). In the QGP the final state phase space of the light particle is 
altered due to the Bose-enhancement (1 + /|j(fc, T)) and Pauli blocking (1 — fp{k, T)) 
factors for gluons and quarks, respectively, recall Eqs. (I2.14p . (I2.15p and (I2.16p . 

B. Charmonium Dissociation in the QGP 

Let us start with the traditional gluo-dissociation process (\E' + (? — > c + c) as illus- 
trated in Fig. [iHl proposed by Bhanot and Peskin in the 1970s [52]. They evaluated 
the interaction between a charmonium state and a gluon within an Operator Prod- 



51 




9 



Fig. 18. Diagrams of the gluo-dissociation process. A charmonium absorbs a gluon 
and is subsequently dissociated. 

uct Expansion (OPE) formalism. The leading order in as turns out to be a dipole 
interaction (r ■ E) between the \l/ and the gluon. For J/ip dissociation cross section 
it results in the following expression, 



-.J/AE,) = f{j) (^) 



ml {Ejesf 



(S.t 



where the rric is the mass of the charm quark, Eg is the energy of incoming gluon 
and tB is the binding energy of J/V", see Eq. (13. ip . We note that the cross section 
agji^{Eg) exhibits a pronounced peak at Eg = ye^, see Fig. [191 Note that since we 
adopt a massive gluon with mg\T=3ooMeV — 580MeV (see Eq. (13.71) ). the peak appears 
at gluon momentum pg ~700MeV for vacuum J/ip binding energy of 640MeV. The 
limitations of the gluo-dissociation process include: (1) As mentioned in Ref. [52] the 
OPE procedure is valid only if the energy scale of the incoming gluon is much less 
the charmonium binding Eg <^ eb- (2) Since Eq. (13.81) is only to the leading order in 
as, it does not include the interaction between J /ip and light quarks. 

Indeed, with in-medium binding energy eb <T the gluo-dissociation mechanism 
turns out to be not numerically efficient in destroying J /ip. To see this we calculate 
the dissociation rate, r,j,, which is a convolution of the density of incoming particles 
(gluons in this context), fg{Pg), with the dissociation cross section, (7g^{s) and the 
relative velocity Vrei between the incoming particle and the charmonium (cf. Eqs. (12. 4p 
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Fig. 19. Comparison of the parton- induced dissociation cross section. Solid line: 
quasifree dissociation mechanism; dashed line: gluo- dissociation with vac- 
uum J/ip binding energy; dotted line: gluo-dissociation cross section with 
in- medium J/ip binding energy of llOMeV; dot-dashed line: gluon thermal 
distribution p2/^(p^) at T=300MeV. 



and dH), 



T^{p,x) 



(Ppg 

(27r)3 



dgfg{Pg)VrelCr 



9*1 



(3.9) 



Here the gluon momentum distribution fg{pg) follows from Eq. (13. Sp . 

Figure [19] shows that the inefficiency of the gluo-dissociation for J/ip with in- 
medium binding energy originates from the peak structure of the gluo-dissociation 
cross section. Note that since we use a massive gluon with mg|T=3ooMeV —580 Me V, 
the gluon energy at Pg=0 is already significantly larger than YeB|T=3ooMeV~160MeV, 
therefore only the "tail" of the gluo-dissociation cross section appears in Fig. [191 
leading to a small overlap between the thermal gluon distribution p^fg{pg) and the 
dissociation cross section, agj/^{pg). The gluo-dissociation cross sections for Xc and 
i/j' show similar peak structures, resulting in small in-medium dissociation rates. 

On the other hand we know that in the limit of incoming gluon energy large 
compared to the charmonium binding energies, e^ > Eg, the interaction between 
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Fig. 20. Diagrams of the quasifree process. A charmonium can be dissociated by either 
a quark or a gluon. 

the c and c inside the charmonium, which occurs on a time scale of l/e^ > l/Eg, 
cannot interfere with the interaction between the incoming gluon and the c{c). In 
this limit the g — "^ scattering cross section should approach the value of 2agc, the 
sum of the probability of the scattering between the gluon and one of the constituent 
charm quarks. Therefore the peak structure of gluo-dissociation cross section will be 
superseded by continuum-hke "quasifree" dissociation cross section. 

Following this idea an alternative (quasifree) dissociation mechanism, i + \1' — )■ 
i + c + c{i = g,q,q), was proposed to describe the inelastic scattering between the 
charmonium and partons [37]. In these processes the incoming parton (gluon or 
quark) collides with the c or c quark in the "$. Since in the QGP charmonia are loosely 
bound states one such collision would be enough to dissociate the \l/, as illustrated in 
Fig. EHl The transition matrix element for the quasifree dissociation can therefore be 
factorized as 



\M^,i^cci{Pi, p^;Pi, Pc, Pc)\'^S''^'' {pi + p^ - Pc - Pc - Pi) 

2 



-^2 



Mci-^ci Pi, — p^;Pi,Pc 



^(4) rn^ 






m* 



P^+Pi-Pc- P, 



X (27r)3 (2Eg) 6^'^^ 



TTlr 



m^ 



-P^ - Pi 



(3.10) 



where | A^^j^ccil is the quasifree dissociation transition matrix element, and \Aic 
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the transition matrix element for the elastic scattering between i and c, calculated in 
Ref. [53]. The momenta pi, p$, pi, pc and pc correspond to the initial state parton i, 
initial \E', final state parton i, final state c and final state c. The overall factor of 2 is 
accounting for the fact that the elastic scattering with either c or the c can dissociate 
the \E'. 

To account for the leading kinematic correction from the residual binding energy, 
the incoming parton needs to be energetic enough to break up the bound state, 
which sets a lower limit for the incoming parton momentum. Overall 4-momentum 
conservation for the process z + ^— T-z + c + cis maintained by assigning the binding 
energy to a decrease in mass of the initial-state charm-quark, c', i.e., rric' = rric — e^. 
In addition, we have introduced a Debye mass, mu = gT, into the denominator of t- 
channel gluon-exchange propagator, 1/t — )■ l/(t — r?7,|,), to regulate the divergence for 
forward scattering (the strong coupling in mo is taken consistently with the coupling 
constant a^ used for the quasifree process). 

Inserting the quasifree transition probability Eq. (13.101) into Eq. (12.151) we can 
evaluate the corresponding dissociation cross section and rate. The comparison of 
the dissociation cross sections between the quasifree dissociation processes and the 
gluo-dissociation process is shown in Fig. [191 In contrast to the gluo-dissociation 
cross section the quasifree dissociation cross section saturates for partons with large 
incoming momentum {k ^ e^). Since the quasifree process may not be the only 
dissociation mechanism for \E', for practical applications we effectively parameterize 
other dissociation mechanism into the quasifree processes by using the strong coupling 
constant a^ (figuring into \Aici-^ci\'^) as an adjustable parameter. We adjust it to 
the J/ip suppression data measured at SPS/RHIC, see Chapter |VT] for the detailed 
procedure. The resulting value of a^ turns out to be ~0.3, quite compatible with 
the short-distance (color- Coulomb) term in the effective potential used to extract the 
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Fig. 21. Temperature and momentum dependence of \1/ dissociation rates. Left panel: 
Temperature dependence of dissociation rates (at p=0) for J/ip and Xc- Solid 
line: J/ip quasifree dissociation rate in the strong binding scenario; dashed 
J/ip quasifree dissociation rate in the weak-binding scenario; dot-dashed 
Xc quasifree dissociation rate in the strong-binding scenario; double- 
dot-dashed line: Xc quasifree dissociation rate in the weak-binding scenario; 
dotted line: J/ip gluo-dissociation rate in the strong binding scenario with 
as=0.3. Right panel: Momentum dependence of J/ip quasifree dissociation 
rates for strong binding scenario (solid line); weak binding scenario (dashed 
line) and gluo-dissociation rate with vacuum J/ip binding energy (dot-dashed 
line) . 



binding energies by the T-matrix potential model. Note that the strong coupling 
constant in the quasifree process is decoupled from that figuring into parton thermal 
masses in Eqs. (13. 5p . (13. 6 p and (13. 7p . 

The temperature dependence of the quasifree dissociation rates for the J/ip and 
Xc in the QGP is plotted for p=0 in the left panel of Fig. [211 ^ind for the J/ip as a 
function of p at selected temperatures in the right panel. 

In the weak-binding scenario, there is rather little difference between the dissoci- 
ation rates of J/tp and Xc especially above T = 200 MeV. Only in the strong-binding 
scenario the larger J/ip binding energy makes a large difference, suppressing its de- 
struction by, e.g., a factor of ~5 at T ~200 MeV relative to the Xc and ip' (not shown); 
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this difference becomes larger (smaller) at smaller (larger) T. For comparison we also 
calculated the rate due to the gluo-dissociation mechanism employing the expression 
derived in Ref. [52] with the same as —0.3 as in the quasifree rate and with e^ ob- 
tained from the strong-binding scenario (note that the Coulombic binding is much 
smaller), which turns out to be inefficient for dissociating J/ip's (and even more so 
for the excited states) and is thus neglected in the following. We also note that to 
achieve a comparable dissociation rate, say, F j/^ ~ 100 MeV at T ~ 300 MeV in a 
hadronic medium, one needs a hadron density of about 5/fm^ for typical thermal 
averaged hadronic dissociation cross section of around Imb [M] • 

The 3-momentum dependence of the rates shows a monotonous increase with 
increasing p, which becomes more pronounced with increasing binding energy (for 
larger e^ a finite 3-momentum facilitates the break-up since, on average, a larger 
center-of-mass energy is available in the collision of the bound state with thermal 
partons). This increase is a simple kinematic consequence of a monotonously increas- 
ing (or even constant) cross section with finite threshold and an increasing parton 
flux encountered by a moving J/ip. 

C. Charm-Quark Spectra and Charmonium Regeneration in the QGP 

In the previous section we discussed in-medium dissociation of charmonium states 
on the microscopic level. The inverse process, the regeneration of charmonium states 
from charm quarks, should also occur and obey the principle of detailed balance with 
the dissociation processes. The expressions for the regeneration rates, Eq. (12.101) and 
Eq. (I2.16p . are similar to those of the dissociation rates, Eq. (12. 4p and Eq. (12.151) . 
The only difference is that the charm-quark phase space distributions, fc(c){x,pc,t), 
appear in the initial phase space. Below we first give a brief review of the typical 
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time evolution of fc(c){x,Pc,t): 

Since the c and c quarks are produced in hard coUisions their initial transverse 
momentum spectra are rather "hard", with (pf) ^3GeV^. (As an estimate for a 
typical temperature of T=250 MeV in heavy-ion collisions, the thermal momentum 
spectra have (p^) ~1 GeV^, much smaller than the initial spectra from hard collision.) 
Later they collide with particles in the (thermalized) medium and gradually equili- 
brate their momentum with the heatbath. If the medium does not cool down they 
will eventually be thermalized and their momentum spectra follow the thermal spec- 
tra, fc{c){pc) oc e~^''''^ with Ec = \/m1 + p^ . However, compared to light quarks the 
charm quark thermalization is a slow process, "delayed" relative to light quarks by a 
factor of ~ nic/T ^ 5. Since the thermalization time for the bulk medium is of order 
~ Ifm/c, the thermal relaxation time for charm quarks is expected be on the order 
of ~ 5fm/c [55], which is comparable to the lifetime of QGP at RHIC. Therefore, 
during the lifetime of the hot medium the charm-quark spectra are expected to be 
partially thermalized. Quantitative calculations of charm quark thermalization can 
be performed with a Boltzmann transport equation for charm quarks. 

The Boltzmann equation for charm quarks is an intego-differential equation, 
which is numerically difficult to handle. Therefore Fokker-Planck equations and 
Langevin simulation techniques are often employed as approximations to the Boltz- 
mann transport equation to describe the time evolution of charm-quark phase space 
distribution function, see Ref. [55] for a recent review. 

In this work we compare charmonium regeneration rates resulting from two lim- 
iting cases of charm-quark spectra: 
1) fully thermalized charm-quark spectra. In the local rest frame, they have the 
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following form, 



(Pp 



oc e 



-^ml+p-^/T _ ^3_^^^) 



2) spectra from initial hard production (pQCD spectra). We employ a parameteriza- 
tion given in Ref. |56] . 

with y4=0.5GeV, 5=6.8 GeV, a=21. This parameterization is based on the charm- 
quark pt spectra generated in 200 GeV proton-proton (p+p) collisions by PYTHIA [57] . 

3) pQCD spectra (13.121) in transverse direction, and a thermal spectrum in longitu- 
dinal {pz) direction, 

'YT oc ^^^^^^ X e-V^/- , (3.13) 

(Pptdpz N{pz){1 + ^pI+pI/BY 

where 

The thermal distribution in longitudinal direction is to mimic the longitudinal smear- 
ing of the center of the mass momentum due to various momentum fractions, x, carried 
by the two colliding primordial partons in cc production process g{q) + g{q) -^ c + c. 
We denote these spectra as "pQCD+thermal" spectra. 

The three types of charm-quark distribution are compared in the left panel of 
Fig. [22] in terms of the pz-integrated pt spectra. They are all normalized to the 
same total charm pair number. The thermal spectra are the softest {{p'f)=l.l GeV^), 
with most of the yield concentrated at low pt. The 3-dimensionally isotropic pQCD 
spectra ((]9^)=3.5 GeV^) and the transversely pQCD + longitudinally thermal spectra 
((Pt)=3.0GeV^) are almost of comparable hardness in the transverse plane, with the 
2+1-dimensional pQCD spectra being slightly softer due to their thermalization in z 
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Fig. 22. Charm-quark pt spectra and J/ip regeneration rates. Left panel: Compar- 
ison of different input charm-quark pt spectra. Right paneh their result- 
ing charmonium regeneration rates evaluated in the strong binding scenario. 
Dashed line: 3-dimensional isotropic pQCD charm-quark spectra; dot-dashed 
line:transversely pQCD + longitudinally thermal spectra; double-dot-dashed 
line: thermal charm-quark spectra with ?7ic=1.73GeV; solid line: transversely 
pQCD with angular correlation + longitudinally thermal spectra. 

direction. 

As discussed in Section UTTllB] the quasifree process is the dominant dissociation 
mechanism for in- medium "$. According to the principle of detailed balance the in- 
verse quasifree process, i + c+c -^ 1 + "^ {i = g,q,q), is the corresponding regeneration 
processes. Its transition matrix element can be obtained by applying detailed balance, 
Eq. f l2.12p . to the quasifree dissociation matrix element, Eq. fl3.10p . 



\Mcci^ii,i{pi, Pc, Pc] Pi, PM/) P'^^'''' {Pc +Pc+Pi-Pi 



p*, 



->2 



dcdc 



Mci^ci Pi,Pc; Pi, — p* 



rrir' 



6^^'^ [Pi+Pc-Pi -p^ 



X 



{2^f{2E,)5^^(p,-^p^) . 
\ m^ J 



(3.15) 



Here, \Mc 



is the transition matrix for the elastic scattering between i and c. 
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The momenta of initial state parton i, initial state c, initial state c, final state i and 
final state \1' are pi, Pc, Ps, Pi, and p^,. rric' = nic — €3 is the mass of the charm quark 
in the final state of the elastic scattering process. Note that the factor of j^ comes 
from summing over final state color-spin degeneracy and averaging over initial state 
degeneracy. In the quasifree approximation the 3— )'2 coalescence process is reduced 
to a 2— 7-2 (quasi-elastic) scattering process. 

We proceed to calculate the charmonium regeneration rate by plugging the 
charm-quark spectra and the transition matrix element, Eq. fl3.15p . into Eq. fl2.16p . 
The resulting J/ip regeneration rates, I3^,{pt)., from different charm-quark spectra are 
compared in the right panel of Fig. [221 The (3qi{pt) can be interpreted as the regen- 
erated number of \E' with a given pt per unit volume per unit time. Here we have 
integrated over the p^ dependence of (3q,{p). 

As expected, the pt dependence of the regeneration rate from thermal (pQCD) 
charm spectra follows a thermal (pQCD) trend. The thermal charm spectra, with the 
largest charm-quark phase overlap in the low-pi region, lead to the largest inclusive re- 
generation rate, f3tot = f d?P P{j>)- The 3-dimensional isotropic pQCD spectra and the 
transverse pQCD + longitudinal thermal spectra lead to significantly smaller inclu- 
sive regeneration rates, amounting to 28% and 47% of that from thermal charm-quark 
spectra, respectively. This is similar to what is found in Ref. |SB], where the dynamics 
of cc coalescence was encoded in a Gaussian Wigner function and the resulting number 
of J /ip from pQCD charm-quark spectra is smaller than that from thermal spectra by 
a factor of 3. However in Ref. [lO] where the (inverse) gluo-dissociation mechanism 
was employed for J /^> regeneration it is found that the inclusive yield of regenerated 
J /%l) from pQCD charm spectra is quite comparable with that from thermal spectra 
(within ~30%). Further investigations are needed to clarify the discrepancy. 

For pQCD-|-thermal charm spectra we also consider a possible angular correlation 
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between c and c in momentum space: it is expected from pQCD that in initial hard 
colhsions back-to-back charm pair production is favored. We study its consequence on 
J/V' regeneration by including a schematic ansatz oidNcc/dO^ (1 — cos^) in Eq. f l2.15p 
(with 6 being the relative angle between the pt of c and c in the transverse plane). It 
turns out that the angular correlation of c and c significantly reduces the regeneration 
rate for high pt J/ip as seen in Fig. [221 The reason is that, compared to low pt J/ips, 
high Pt ones are more likely to be regenerated from a cc pair with small angle between 
them, the probability of which is suppressed by the angular correlation factor (1 — 
cos 6*). The inclusive regeneration rate from the angular correlated pQCD+thermal 
spectra is 11% of that from thermal charm spectra. We subsequently plug /3*(pt) into 
the Boltzmann transport equation and the numerical results for the inclusive yield 
and Pt spectra of the regenerated J ftp will be compared in Chapter I VII 

D. Charmonium Dissociation in Hadronic Matter 

As the fireball expands, the medium keeps cooling. When the medium temperature 
decreases to the critical temperature, Tc, the deconfined QGP undergoes a phase 
transition to the confined hadronic phase. The charmonium dissociation rate in the 
hadronic phase is expected to be small compared to QGP due to the smaller light par- 
ticle density and the larger charmonium binding energy. However, for a quantitative 
calculation of charmonium yields in heavy-ion collisions, dissociation/regeneration in 
hadronic phase needs to be taken into account. 

The microscopic approaches for the hadronic charmonium dissociation can be 
divided into two categories, based on either quark or hadronic degrees of freedom. 
The approaches based on hadronic degrees of freedom are often less cumbersome yet 
very effective in assessing many dissociation processes, see Ref. [59] for a recent review. 
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Fig. 23. J ftp dissociation rates by pions and rhos as a function of the temperature T 
of the hadron gas. Figure taken from Ref. [61]. The full line is the sum of the 
three contributions, J/ip + v" (dashed line), J/ip + p ^>- D* + D* (dot-dashed 
line) and J/ip + p ^>- D + D (dotted line). 

In this work we adopt dissociation rates based on an effective meson Lagrangian with 
a local SU{A) flavor symmetry [5H|60l[6T]. The considered interactions are J/ip with 
pions and rho mesons which are the most abundant particles in the hadronic medium. 



p+J/i) ^D + D 
p + J/i/j -^D* + D* . 



(3.16) 
(3.17) 
(3.18) 



We neglect the interactions between J/ip and kaons (In Ref. [62], the total cross 
section of J ftp dissociation by kaons is found to be much smaller than that by pions). 
The resulting dissociation rates are shown in Fig. [231 with the main contribution 
given by the p + J/ip — )■ D* + D* process. The vacuum masses are assumed for all 
the hadrons including J/ip. 

For dissociation rates of Xc and ip' we assume a geometrical scaling |114] of the 
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dissociation rates for J /ip by the respective charmonium radii, namely, 

r?^^ = (;^) ' r J4^ ^ 2.36 r J// , (3.19) 

T^/^ = (^) ' r^^/ ^ 3.73 rf// . (3.20) 

The radii of the excited charmonia states are estimated by non-relativistic potential 
models from Ref. [24] . 

As expected, the overall charmonium dissociation rates in the hadronic phase are 
significantly smaller than in the QGP phase, resulting in a rather mild charmonium 
suppression in the hadronic medium compared to the QGP. As in the QGP, \1' can also 
be regenerated in the hadronic matter [63]; their regeneration rates can be obtained 
from their dissociation rates using the principle of detailed balance. 

E. Charmonium Spectral Function 

So far we have obtained the charmonium binding energy e^ and in-medium charm 
quark mass m* (from the potential model) and its dissociation rate (from quasifree 
approximation in QGP and S'f/(4) effective theory in HG). The charmonium binding 
energy e^ and in-medium charm quark mass ?7i* allow us to infer the charmonium 
pole mass via m^=2m*-eB. These pieces of information figure into the charmonium 
spectral function, Eq. (I1.20p . which is the imaginary part of the Fourier transform 
of the charmonium current-current correlation function, Eq. (ll.lOp . Due to current 
limitations mentioned in Section UllDl we can not reliably extract the charmonium 
spectral function from the charmonium current- current correlation function calculated 
by lattice QCD. However, we can still utilize the rather precise correlator ratios, 
Eq. f ll.24p . as a constraint on the consistency among the charmonium pole mass, 
binding energy and its dissociation rate. 
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To be specific we adopt the following strategy: We "reconstruct" in-medium 
charmonium spectral functions using a relativistic Breit-Wigner + continuum ansatz, 
where the "^ width and mass figure into the Breit-Wigner part while the continuum 
is determined by the open-charm threshold (2m*). For a more realistic evaluation, 
we include a polestrength factor, Z^(T), for the Breit-Wigner strength and a non- 
perturbative rescattering enhancement in the continuum |6H ES]- The vanishing of 
the polestrength factor furthermore serves to estimate the dissociation temperature 
of the ground state in each channel. 

We first construct a model spectral function in vacuum, consisting of a zero-width 
bound-state and a perturbative (leading order) continuum part. 



(T*(w) = A^ S{uj - m$) H —^Q{uj - yi^)cj^wl ^(a + b-^) . (3.21) 

Here, Nc=3 is the number of colors and the coefficients (a, b) = (1, —1), (2, 1) charac- 
terize the scalar and vector channel, respectively [66]. The open-charm threshold in 
vacuum, ^/so, is assumed to be given by twice the free D-meson mass, ^/so = 2mD = 
3.74 GeV. The coefficient A^^, is related to the overlap of the wave-function, /?j/^(0), 
or its derivative, i?^^(0), at the origin [66| [67]. 

Aj/, = ^-^\Rj/m? , \. = ^\KM\' ■ (3-22) 

Xc 

These quantities can be estimated from the electromagnetic decays widths via [67] 

where a=l/137 is the electromagnetic coupling constant and eg = 2/3 the charge of 
the charm quark (we use ree=5.55keV for the J/ip and r^^=2.40keV for the Xco)- 
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Fig. 24. Temperature dependence of the strength of the resonance part of the S-wave 
spectral function, Z^{T), for the strong (sohd hne) and weak (dashed hne) 
binding scenarios. 



The resulting relations between Aij, and r^_ 
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(3.24) 



The free J/tp and Xco niasses are taken at their empirical vacuum values. The coef- 
ficient B^j, in the continuum part of Eq. (I3.2ip equals to one in the non-interacting 
limit. To account for rescattering, which is particularly important close to threshold, 
we scale it up to match the continuum as calculated from the vacuum T-matrix in 
Ref. [50], amounting to -Bj/ij, ^ 2 and B^^ :i^ 4 in the vector and scalar channel, 
respectively. For simplicity we neglect ip', x'c ^^^1 higher excited states which play 
little role in the correlator ratios, Eq. (I1.24p . 

At finite temperature we replace the (5-function bound-state part by a relativistic 
Breit-Wigner (RBW) distribution while the continuum part is assumed to be of the 
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Fig. 25. Spectral functions in the vector channel for the strong (left panel) and weak 
(right panel) binding scenarios. 



same form as in the vacuum, 



aii,{u) = Ail, Z^i{T) 



2w 



wr*(T) 



vr {uj^-ml{T)Y + uj^Ti,{Ty 



+ 



B<i,N, 
87r2 



■Q{u - V^))u^Jl - ^(a + 6^) . 



(3.25) 



The in-medium continuum edge, s(T), is now taken as the charm quark threshold 



at finite temperature, \/s{T) = 2m* (T), consistent with the potential model, see 
Fig. [161 The RBW term includes: (i) the in-medium charmonium mass, m^lT), 
extracted from Eq. (13. ip based on Figs. [16] and [T71 (ii) the width F^ identified with 
the inelastic dissociation width discussed in the previous section; (iii) the aforemen- 
tioned polestrength factor, Zq,(T), representing the modification of the strength of 
the bound-state part at finite temperatures relative to its vacuum value (v4,j,), with 
Z^\t=o=1. The Z<i,{T) is adjusted to minimize the deviation of the correlator ratios 
from one. 

The resulting 2',j,(T) for J/ip (vector channel) is plotted in Fig. [24l from which we 
extract its dissociation temperature Tjf^=2.0{1.25)Tc in the strong (weak) binding 
scenario. Similar analysis in the scalar channel yields Xc dissociation temperatures 
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Fig. 26. Ratio of vector channel correlator to the reconstructed correlator for the strong 
(left panel) and weak (right panel) binding scenarios. 
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Fig. 27. Spectral functions in the scalar channel. The left (right) panel is for the strong 
(weak) binding scenario. In the weak-binding scenario Xco has already melted 
at Tr. 
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of T°^^^=1.3(1.0)Tc in the strong (weak) binding scenarios. We assume that Xci and 
Xc2 have the same dissociation temperatures as the Xco- For ip' we simply assume its 
dissociation temperature to be Tc for both the strong- and weak-binding scenarios. 

To comprehensively illustrate the medium effects we plot the final spectral func- 
tions for the vector channel in the strong- and weak-binding scenario in the QGP in 
Fig. [251 and their corresponding correlator ratios in Fig. [221 the spectral functions 
for the scalar channel are displayed in Fig. [271 

We see that the correlator ratios are indeed close to one, as found in IQCD [311 
[301 [29] . In the hadronic phase (not shown), we assume vacuum masses for both 
charmonia and open-charm hadrons, which automatically ensures that the correlator 
ratios are close to one (deviations due to small charmonium widths in hadronic matter 
are negligible). 
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CHAPTER IV 

THERMAL FIREBALL DESCRIPTION OF MEDIUM EVOLUTION 

To solve the kinetic equation for the charmonium {^=J /ip, Xci 4^') evolution we 
need a description of the expanding medium in heavy-ion collisions. In particular, 
the temperature and the volume of the medium at any given time are required for 
estimating parton density and charm quark density, respectively, which are essential 
components for calculating the \Ef dissociation and regeneration rates. The traditional 
approaches describing the medium evolution include transport models, see Refs. [68l 
EH] for reviews, and relativistic hydrodynamics, see [70] for a review. In this work we 
however describe the medium evolution with a fireball model [7T| [72] which is simpler 
yet captures the basic features of relativistic hydrodynamics. In Section lIVIlAl we 
discuss the fireball description of the spatial expansion of the medium. In Section HVllB] 
we discuss the equation of the state of the medium and extract the temperature of 
the fireball. In Section IIVO we review the popular blastwave formula for estimating 
particle pt spectra with local thermal distributions boosted by the collective flow of 
the expanding source. 

A. Fireball Expansion Profile 

The fireball model approximates the medium created in heavy-ion collisions as a 
boost-invariant thermal fireball. The "boost-invariance" originates from the exper- 
imental fact that the rapidity (y) distribution of the charged particles, dNch/dy, is 
constant in the mid-rapidity {y ~ 0) region. This means that the medium in the cen- 
tral region is invariant under Lorentz boosts in the longitudinal {z) direction, which 
further implies that all thermodynamic quantities characterizing the central region 



depend only on the longitudinal proper time T=y/t'^ — z'^, recall Fig. [71 
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The fireball volume, Vpb(t), expands cylindrically according to 

VFB{r) = {zo + v,r + -a^) tt tq + ^^ ^ 4.1 

2 y a± J 

with expansion parameters Vz,az, a-± chosen so that the results are consistent with the 
experimental data on the final light-hadron fiow and resemble the evolution by hydro- 
dynamical calculations. The relativistic form of transverse acceleration (a/1 + a']_T'^~ 
l)/a± limits the surface speed, Vs{t), to below the speed of light in the large r limit. 
For small r it recovers the non-relativistic form with constant transverse acceleration, 
'Vs{t) ~ a±T. Note that the fireball expansion profile, Eq. (14. ip . depends only on the 
longitudinal proper time, T=\/t'^ — z'^, and not separately on t or z, as required by 
the boost invariance of the medium. 

The initial transverse radius tq represents the initial transverse overlap of the two 
colliding nuclei at a given impact parameter b, while the initial longitudinal length, 
Zq, is related to thermalization time Tq through Zq ~ Ayro where Ay=1.8 represents 
the typical longitudinal rapidity coverage of a thermal fireball. We assume that at a 
formation time of ro=1.0 (0.6) fm/c the medium at SPS (RHIC) first thermalizes with 
all the entropy, Stot{b), being built up. The latter is estimated from the multiplici- 
ties of observed charged particles and assumed to be conserved during the adiabatic 
expansion. 

B. Equation of State of the Medium 

To determine the temperature of the system we utilize the equation of state (EoS) 
of the medium. An important insight from the success of the ideal hydrodynamic 
description of the medium evolution is that the total entropy of the system is con- 
served. Together with the information of the fireball expansion profile the entropy 
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Fig. 28. Comparison of energy density from an ideal massive parton gas with data 
from lattice QCD [S]. Dashed (dot-dashed) line: energy density (pressure) 
calculated from an ideal massive quark and gluon gas with (7^=2.3. 

density of the fireball at each given moment can be obtained. This quantity allows 
us to estimate the temperature of the fireball through comparison with the thermo- 
dynamic equation of state (the entropy density as a function of the temperature, T). 
In QGP we model the medium by an ideal gas of massive quarks and gluons, with 
the entropy density as a function of temperature given by 

where di is the color-spin degeneracy of partons, the fi{k;T) are the thermal distri- 
butions for massive quarks and gluons, given by Eqs. (13. 4p and (13. 3p respectively, and 
the plus sign in "±" is taken for gluons {i=g) and the minus sign is taken for quarks 
{i=q,q)- For the temperature dependence of thermal quasiparticle masses irtiiT) we 
take guidance from pQCD calculation, see Eqs. (13. 5p . (13. 6p and (13. 7p . The strong cou- 
pling constant g is estimated by matching the resulting energy density of the parton 
gas to IQCD calculations, see Fig. |28l resulting in Qs ~2.3. Note that this g is decou- 
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RHIC-200 mid-y: s/pg=250 
SPS-17.3: s/pg=26 
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Time (fm/c) 



Fig. 29. Time profiles of temperature for central collisions of heavy nuclei (participant 
number A^part=380) at RHIC (^5=200 AGeV; solid line: mid rapidity; dashed 
line: forward rapidity) and SPS {y/s=17 .3 AGeV; dot-dashed line). 

pled from the a^ in the quasifree dissociation cross section. In the hadronic medium 
an EoS similar to Eq. (14. 2 p is employed but including 76 mesonic and baryonic states 
up to masses of 2 GeV. For particles which do not decay strongly, e.g., pion, kaon, 
antiproton, we include their corresponding chemical potentials to maintain their abun- 
dances when the system cools down from chemical to thermal freeze-out [7T]. The 
critical temperature, Tc=170 (180) MeV at SPS (RHIC), is roughly consistent with 
thermal- model fits to observed particle ratios [73] and predictions of lattice QCD |74] . 
A freeze-out temperature of Tfo ~ 120 MeV terminates the evolution and results in a 
total fireball lifetime of rfo=10-12fm/c for central A+A collisions. The resulting tem- 
perature evolution as a function of time r is displayed in Fig. [29] for SPS and RHIC. 
Note that there is little difference between mid {\y\ < 0.35) and forward rapidity 
{\y\ G [1.2,2.2]) for Au+Au collisions at RHIC due to the slowly varying rapidity 
density of charged particles over this y range [75] , cf . Ref . [76] for more details. 

For heavy-ion collisions with lower center-of-mass (cms) energies, e.g., at the 
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Fig. 30. Excitation function of the baryon chemical potential ^b{Tc) at the assumed 
phase transition line. The points correspond to values extracted from particle 
ratios. 

future FAIR facility, stronger stopping of incoming nuclei entails the medium to be 
more asymmetric in terms of the net baryon number, corresponding to a larger baryon 
chemical potential ^Bi see Fig. |30l and lower initial temperature, T. In this regime the 
thermodynamic properties of the medium are characterized by both T and /^b(^ T). 
In order to determine the temperature evolution in such systems we resort to a second 
conserved quantity during the fireball evolution, which is the net baryon number, Nb- 
With both baryon density, ub, and entropy density, s, given at each time we are able 
to solve for ^b{t) and T{t) by using 



fi{k;T,fXB) log fi{k;T,fXB)] 



nB{T,^B) = Y^ 



di 



f 3 X (27r)3 



d'k[{f,{k-T,^B)-f^{k-,T,^B)] 



(4.3) 
(4.4) 
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Fig. 31. Temperature evolution of heavy-ion collisions at FAIR. Left panel: Trajectory 
of heavy-ion collisions on the QCD phase diagram at FAIR (y^=8.8AGeV). 
Right panel: Time profile of temperature for central collisions of heavy nuclei 
(participant number Npart =380) at FAIR {y/s=8.8AGeV). 



The thermal quark distributions at finite fj,B are 



fq{k]T,fiB) 



fq{k- T, hb) 



1 



exp fV^I^ztE/1 ) + 1 



exp 



/^/^+^+^^\ _^^ ■ 



(4.5) 









As part of the EoS of the system, the phase boundary fisiTc) (instead of a single Tc) 
is determined with guidance from thermal-model fits to observed particle ratios at 
various center-of-mass energies [73], see Figs. [30] and [3T1 

The resulting time evolution of the fireball is characterized by a trajectory in the 
QCD phase diagram at fixed s/ub ratio, see the left panel of Fig. EH We note the 
zigzag structure of the trajectory, which reflects the fact that the medium is reheated 
by the latent heat when it undergoes the phase transition, as first pointed out in 
Ref. [77]. 

The resulting temperature evolution as a function of time r for central Pb+Pb 
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collisions at FAIR energies is displayed in the right panel of Fig. [HH The thermaliza- 
tion time tq is assumed to be 1.3fm/c at FAIR energies. With an initial temperature 
lower than 160 MeV the QGP phase still lasts for l-2fm/c, which is followed by the 
mixed phase with duration of around 2fm/c. The overall fireball lifetime is again 
around lOfm/c. 

C Blastwave Description of Charmonium pt Spectra 

In addition to the temperature evolution profile the thermal fireball model allows to 
define a flow field, which enables us to estimate the transverse momentum [pt) spectra 
of locally (kinetically) thermalized particles, boosted by the flow field. Historically 
this analysis was developed by Siemens and Rasmussen [7S], and is referred to as 
the blastwave model. Later the blastwave model in a boost-invariant medium was 
developed by Schnedermann, Sollfrank and Heinz [79]. Here we follow Ref. [79] and 
briefly review the derivation of the "blastwave formula" . 

Since the medium is not static in the lab frame we employ the Cooper-Frye 
formula [HD] to count the particle number in the freeze-out hyper-surface, 

^1^ = (^ /^Sm(-K/(-,p) = ^. j d^^^-)P'^-''"" ■ (4-6) 

Here the S^ is the freeze-out hyper-surface, u'^ is the four velocity of a given medium 
cell, and f{x^p)=de~^'^/^ is the thermal particle distribution in the local rest frame 
with d being the degeneracy factor. 

Since the freeze-out is commonly assumed to occur at fixed longitudinal proper 



time, T=\Jt'^ — z^, it is convenient to parameterize the freeze-out hyper-surface by 

x^ = {tjo cosh r], r cos (j), r sin 0, tjq sinh rj) , (4.7) 
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where r]=tanh ^{z/t) is the space-time rapidity. The volume element of the hyper- 
surface is therefore 

dT,^ = Tfo {coshrj, 0, 0, sinhr^) r dr dr] dcj) . (4.8) 

The particle four-momentum can be parameterized as 

p'^ = {rrit cosh y, pt cos 0p, pt sin 0p, rrit sinh y) , (4.9) 



where |/=tanh [pz/ E) is the momentum rapidity of the particle and mt=\/m? + p^ 
its transverse mass. Therefore we have 

p ■ dT, = TfoiTit cosh(?/ — ri)rdrd'r]d(f) . (4-10) 

A parameterization of the four-velocity m^ of a medium cell can be written as 

u'^ = (cosh p cosh y^*^ , sinh p cos (/)p, sinh p sin 0p, cosh p sinh y'^^ ) . (4-11) 

Here y^^'"'' = tanh~^(f^'^'') is the longitudinal rapidity of a medium cell. According 
to boost-invariant condition (f^^" = z/t), it is equal to the longitudinal space-time 
rapidity, namely, y'^^" = 77; p=tanh~ (t>^^") is the transverse rapidity of the medium 
cell. Therefore the four product p ■ u can be expressed as 

p ■ u = miCoshpcosh(r7 — y) — Pi sinhpcos(0 — </)p) . (4-12) 

Plugging Eqs. (14.101) and (14.121) into Eq. (14.61) we obtain the one particle momentum 
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spectrum in the lab frame 



^~KZ ~ ToZVa / "r / dr] rdr 
'o J~oD Jo 



'd^p (27r;^3 



X TfoT^t cosh(y — rj) exp[(—mt cosh p cosh (r] — y) + Pt sinhpcos(<^ — (f)p))/T] . 

(4.13) 

After the integration over rj and (j) is performed we arrive at the blastwave formula, 

^dN dN d f J r^ / ""^i cosh p Pj sinh p ^ 



E- 



f , :r^ .mtcoshp. ^ .pjsinhp, ,.-... 

;Tf,m,J rdrK,{ ' ^ ^ )/,(^l__J^) . (4.I4) 



d'^p dyd^pt Iti'^ 

Here K\ and Jq are modified Bessel functions. The transverse rapidity of a medium 
cell p=tanh~^(f^^") is given by the fireball model by assuming a linear transverse flow 
profile with 

<\^ = ^v, , (4.15) 

where R is the radius of the fireball and Vs{T)=a±T is the transverse flow velocity at 
the surface of the fireball. 

In Chapter I VI I we will utilize the blastwave formula to estimate the transverse 
momentum spectra of \l/ regenerated from thermal charm quark spectra. 
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CHAPTER V 

PRIMORDIAL CHARMONIUM PRODUCTION 
The central purpose of studying charmonium production in heavy-ion collisions is to 
utilize charmonium as a probe of the hot and dense medium. Therefore charmonium 
production in the pre-equilibrium stage is necessary as a baseline for assessing any 
modifications due to the hot medium. Since charm-anticharm production is a hard 
process at SPS and RHIC (the production time of charm quark pair, t'^'^ ~0.07fm/c, is 
shorter than nucleus passage time, Tpass ~0.13fm/c at RHIC), it can be approximated 
as superposition of production in elementary nucleon-nucleon collisions. Thus, in 
a first approximation, the initially produced charmonia in A+A collisions can be 
estimated from their production in p+p collisions (scaled by the number of binary 
nucleon-nucleon collisions in an A+A collision). The main purpose of this chapter is 
to apply further corrections specific for p+A and A+A collisions. These corrections 
are usually referred to as cold nuclear matter (CNM) effects. We first give an overview 
of charmonium production in p+p collisions in Section IVlXI Then we proceed to the 
CNM effects in Section |V]|Bl 

A. Charmonium Production in p+p Collisions 

In p+p collisions charmonia are produced in two steps: 

1) The cc pairs are produced through hard collisions (with large momentum transfer) 
between the partons from the two colliding protons. The leading order processes in 
pQCD include the gluon fusion (dominating at high energies) and quark annihilation, 

g + g-^c + c, g + g— )-c + c. (5.1) 
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As hard processes the cc production cross section in p+p colhsions can be factorized 
into a convolution of incoming parton distribution function, /i(x, Q^), and the parton 
scattering cross section, o"'-'^['^^](xi,a;2, Q^), to yield, 

^7Viv^N^g2) ^ l\^^ ['dx2y2fi{x„Q')fj{x„Q')a'^^^''\x„X2,Q') . (5.2) 
Jo Jo ij 

Here, Xi and X2 are the momentum fractions of the two colliding partons within the 
two colliding protons. The large virtuality, Q^, allows us to compute the partonic 
cross-section, a{xi,X2, Q'^), as a perturbative expansion in powers of as{Q^). Accord- 
ing to the uncertainty principle the time scale of the perturbative partonic process is 

~i/g. 

2) c and c quarks which are close to each other in the phase space ( "pre-resonance" cc 
states) may develop into a charmonium state through non-perturbative "final-state" 
interactions. Compared to charm quark production this step is much slower, with a 
typical formation time for charmonium states on the order of 1/e^ ^ l/(2mc). Calcu- 
lations based on a non- relativist ic Schrodinger equation suggest the typical formation 
time of charmonium states to be Tq, ~l-2 fm/c [20]. Due to the large separation in the 
time scale the quantum interference between charm quark production and charmo- 
nium production is suppressed. Therefore the factorization between the charm quark 
production and charmonium production is usually assumed, 

(^A+B^^+X ~ 2_^ / '^'^cc CrA+B^cc[n.]+x($cc, ""^c) Fcc[n]^'^ (^ cc) , (5.3) 

n J 

with a sum over possible cc[n] states and an integration over available cc phase space 
d^cc] F represents a non-perturbative transition probability for a pair of off-shell 
cc to a charmonium state, \E'. The microscopic mechanisms for this transition are 
still under debate. Three widely discussed models in the literature are the color 
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evaporation model (CEM), the color singlet model (CSM) and the non-relativistic 
QCD (NRQCD) model. 

The CEM [8T1 [821 [83] assumes that all cc pairs with invariant mass less than the 
threshold of producing a pair of open-charm mesons, regardless of their color, spin, 
and invariant mass, have the same probability to become a charmonium. That is, 
the -Fcc[„]-5>* in Eq. f l5.3p is a constant for a given quarkonium state, which is usually 
obtained by a fit to the data. 

The CSM [m Uni [861 Ell EH] assumes that only a color-singlet charm quark pair 
with the right quantum number can become a charmonium of the same quantum 
number and the transition from the pair to a meson is given by the charmonium wave 
function overlap at the origin. 

The NRQCD model [671 EH [90] allows every cc[n] state to become a bound 
charmonium, while the probability is determined by corresponding non-perturbative 
matrix elements Fcc[n]-^<ii- The Fcc[n]^^! encoding the non-perturbative dynamics is 
expanded in terms of local matrix elements in a power series of the heavy quark 
velocity, v. 

Each of these models has its advantages and limitations in explaining available 
experimental data, see Ref. [91] for a review. The ^ production mechanisms in 
p+p collisions also have infiuence on \& production in A+A collisions. For example, 
recent RHIC measurements [92l [93] found a reduced suppression for high pt J/ip^s 
compared to low pt ones in A+A collisions, which lends support to the picture that 
the "pre-resonance" cc states are in a color-singlet state, see Section IVI|[B][2][b] for a 
more detailed discussion. 

Since our purpose is to utilize charmonium as a probe to the hot medium we rely 
on experimental values on J/tp and cc production in p+p collisions as the baseline for 
assessing in-medium physics in A+A collisions. For the p+p charmonium production 
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cross per unit rapidity we take the values da pp / dy =37 mh [51] for -\/i=17.3 AGeV 
Pb-Pb [94J (with ca. 40% uncertainty) and da*/rfy=750(500) nb for yi=200AGeV 
Au-Au [95j at mid and forward rapidity (with ca. 10(20)% uncertainty). The input 
charm quark N+N cross section at SPS energy, dcrgc/dy (y=0)=2.2 /ib, is taken 
from a recent compilation of experimental data in Ref. [96] • For full RHIC energy, 
we use dcTgc/dy (y=0)=123±40/ib, in line with recent PHENIX measurements [97] . 
We reduce the input charm quark cross section at forward rapidity by 1/3, dcTpp/dy 
(l/=1.7)=|d(jpp/dy(y=0), according to recent experimental data [98] . 

B. Charmonium Production in A+A Collisions 

After having fixed charmonium production in p+p collisions the next step is to scale 
it by the number of binary nucleon-nucleon (N+N) collisions, A^cou, to estimate char- 
monium primordial production in A+A collisions, which is a standard procedure for 
a hard probe. A^cou is usually calculated with the Glauber model, which plays an 
important role in connecting p+p and A+A collisions. Below we briefly review the 
main results from the (optical) Glauber model which are needed in this work. The 
emphasis is placed on concepts relevant for \1' production, for a more comprehensive 
review of the Glauber model, see Ref. |99l . 



1. Brief Review of Glauber Model 

In this section we give a brief review of the Glauber model and introduce two im- 
portant quantities linking A+A collisions with elementary nucleon-nucleon (N+N) 
collisions: 1) the number of binary collisions, Ncoii, 2) the number of participants or 
wounded nucleons, A^part, which are nucleons from the projectile or the target which 
suffer at least one inelastic collision. In the Glauber model the collision between two 
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Projectile B 



Target A 




"t-A 



Kt 



\j 



a) Side View 




b) Beam-line View 



Fig. 32. Schematic representation of the geometry of Glauber model. Left panel: 
transverse view. Right panel: longitudinal view. The picture is taken from 
Ref. 199]. 



nuclei, A and B, consisting of A and B nucleons, respectively, is considered as a su- 
perposition of (binary) collisions of the individual incoming nucleons. The geometry 
of the Glauber model is schematically sketched in Fig. [321 

The inputs of the Glauber model are 1) nuclear charge densities, usually taken 
as a Wood- Saxon density profile. 



p{r) 



Po 



(5.4) 



1 + exp (^) ' 

where po=0-17fm~^ is the nucleon density in the center of the nucleus, R is the nuclear 
radius and a is the thickness of the nuclear skin. For ^^^Au, i?=6.38fm, a=0.535fm. 
For ^"^^Pb, R=6.62im., a=0.549fm [100] : 2) Inelastic nucleon-nucleon cross section, 
(J^^i, which can be obtained from experimental measurements, e.g., (7,^^=30(42) mb 
at SPS (RHIC) [n]. 

The nuclear thickness function, 



TA{Xt) 



dz pA{xt,z) 



(5.5) 
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can be interpreted as the probability of finding a nucleon within a unit transverse 
area at Xt of nucleus A. The pA is the nuclear density in Eq. (15 ■4p normalized to 1, 
namely, 

PAd^Xtdz = -7 / pAd^Xtdz = 1 . (5.6) 

Therefore TA{xt)TB{xt — b) is the joint probability of finding a pair of nucleons from 
nuclei A and B, respectively, within the common unit transverse area at Xj. Their 
corresponding number of collisions is given by TA{xt)TB{xt — b)a^i d'^Xf. Integrating 
over the transverse plane we obtain the total number of collisions contributed by this 
pair, 

K^nib) = crZ TAsib) = aZ J Ta{x,)Ts{S, - b) d'x, , (5.7) 

where TAB{b) is called the nuclear overlap function. Since from the two nuclei A and 
B a total number of AB such pairs can be found and each of them contributes an 
equal number of collisions, we obtain the following expression for the total number of 
binary collisions 

N,,,,{b)=ABTAB{b)aZ. (5.8) 

The centrality of heavy-ion collisions is often expressed in terms of the number 
wounded nucleons (participants), A^part- In the Glauber model. Apart can be esti- 
mated as follows: The probability for a given nucleon from nucleus A to be located at 
transverse position xt is TA(xt), and the probability for this nucleon to collide with a 
nucleon from nucleus B (located at (xt — b)) is T^(xt) Tsixt — b) (jZ- The probability 
of not colliding is thus TA^Xt) (1 — Tb^Xi — b) 0";^^)- The probability of not colliding 
with any of the B nucleons from nucleus B is thus TA^Xt) [1 — Tb^Xi — b) cr?^^J^. 
Therefore the probability for the nucleon at Xt suffering at least one collision is 
Ta^Xi) (1 — [1 — Tsi^Xt — b) cTin^i]'^). Integrating over the transverse plane we obtain 
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the probability for a given nucleon in nucleus A suffering at least one collision, 

Pl{b) = [ta (ft) \l-\l-TB{xt- b) aH ""I d'xt . (5.9) 



Since there are A{B) nucleons in nucleus A(B) we obtain the total number of wounded 
nucleons (participants) in A+B collisions at impact parameter b as 

N,^,,ib) = A P±{b) + B PZib) 

= a[ta (xt) |l - [l - Tb {xt - 6) aZy] d^^t 

+ B j TB(xt-b)[l-[l-TA {xt) ^a^} (fxt ■ (5.10) 

2. Cold Nuclear Matter Effects 

The notion that ^ production in A+A collisions can be viewed as superposition of 
independent N+N collisions is only approximately true. The deviation of primordial 
\E' production in A+A from A^coU-scaled p+p collisions is usually attributed to the 
so-called cold nuclear matter (CNM) effects. 

In this section we examine the following three aspects of CNM effects: 1) Nu- 
clear shadowing, 2) Cronin effect 3) Nuclear absorption. These CNM effects can, 
in principle, be estimated from p+A collisions where no hot medium is expected to 
form. 

a. Nuclear Shadowing 

It is a well-established fact that the partonic structure of high-energy nuclei is different 
from the incoherent superposition of the constituent nucleons, see Ref. jlOl] for a 
recent review. This modification is usually parameterized by 

Rf{x,Q')= {V^ni, . ^ = <l.<i.9. (5.11) 
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Fig. 33. Modification of parton distribution function in nucleus. Figure taken from 
Ref. [Tn2] . 

defined as the ratio of the parton distribution function for a nucleon inside the nucleus 
(nPDF), //^(x, Q^), to the corresponding one for a free proton (PDF), /j(x, Q^), 
where x is the longitudinal niomentum fraction of the parton within the nucleon. 
As illustrated in Fig. [33] different names have been assigned to these modifications 
depending on the relevant range of x under consideration: 1) Shadowing for the 
suppression observed at small (x < 0.05). 2) Antishadowing for the enhancement at 
moderate values of 0.05 < a; < 0.3. 3) EMC effect ioi the suppression observed in the 
region 0.3 ^ x < 0.7; and 4) Fermi motion for the enhancement when x — t- 1. Let 
us estimate the impact on charmonium production from the modification of nuclear 
parton distribution function by assuming that the initially produced charm quark 
pair has the same rapidity as the charmonium into which they evolve. Then the 
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momentum fraction of the incoming partons x is 

TTlf 

a^i,2 = , exp {±y) , (5.12) 

\/snn 

with the transverse mass rrit = \/vn?^ + p^ of \1/, and y being its momentum ra- 
pidity. In 1/5=17. 3AGeV Pb+Pb colhsions at SPS the relevant x ~0.2 is in the 
ant i- shadowing region. In y^=200AGeV Au+Au collisions at RHIC, for charmonia 
produced at mid-rapidity, the relevant x for partons from both colliding nuclei is 
around 0.02, close to the transition from the shadowing to the ant i- shadowing region. 
For charmonia produced at forward rapidity, y ~1.7, the x of the parton from the 
forward going (in the same direction with the produced \i/) nucleus is around 0.1 
(anti-shadowing region), while the x of the backward going parton is around 0.003 
(shadowing region). In most parameterizations of nPDF, the shadowing of the back- 
ward parton is generally considered to be stronger than the anti-shadowing of the 
forward going parton, see Fig. [331 So the overall effect from the modification of par- 
ton distribution function is expected to cause suppression of \E' production relative 
to p+p collisions. However, for high "pt ^ produced at forward y at RHIC the x 
of the backward going parton could be shifted to 0.01 (for ]9i=10GeV), where the 
shadowing is small and could well be compensated by the anti-shadowing of the for- 
ward going parton. Therefore it is expected that for high p^ "^ produced at forward 
y the suppression due to nuclear shadowing should be suppressed. The shadowing 
/ anti-shadowing effects are experimentally observed in d+Au collisions at RHIC, 
see Fig. [Ml The enhancement (suppression) of charmonia produced at backward 
(forward) rapidity is due to anti-shadowing (shadowing). 
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Fig. 34. Rapidity dependence of J/ip RdAu in d+Au collisions at i/s=200AGeV mea- 
sured by PHENIX [11)8] . 

b. Cronin Effect 

The Cronin effect refers to an enhancement of hadron production at intermediate and 
high pt in P+A relative to p+p collisions (scaled by A^cou)- This effect is generally 
attributed to multiple soft scatterings of the projectile partons propagating through 
the target nucleus before the hard scattering. From the transverse kicks in the soft 
scatterings the partons acquire additional {pf) and therefore the {pf) of the finally 
produced charmonia increases correspondingly. Fig. |35] illustrates the Cronin effect 
at SPS energies, where the (pi) of charmonium produced in p+A collisions increases 
with size of the colliding nucleus (approximated by the path length travelled by the 
colliding parton). The Cronin effect for J /if) is currently difficult to quantify at RHIC 
energies due to large uncertainties associated with d+Au data. 
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Fig. 35. I\pf) of the J/^ as a function of the geometric length of matter, L, traversed 
by partons in the initial state. The figure is taken from Ref. |104j . 

c. Nuclear Absorption 

In p+A or A+A collisions the "pre-resonance" cc states (those cc pairs close to each 
other in phase space which would form charmonium if there were no rescattering 
off surrounding particles, ie., in p+p collisions) are subject to dissociation through 
inelastic collisions with passing-by nucleons before they are fully developed into char- 
monia. Nuclear absorption is observed in p+A collisions at SPS, where the energy 
deposited is too small to create a hot medium: charmonium production is substan- 
tially suppressed relative to ( A^cou scaled) p+p collisions, see Fig. |36l At RHIC energy 
nuclear absorption is observed in d+Au collisions, see Fig. [STJ which is weaker than 
at SPS, presumably because the time scale over which the two nuclei pass through 
each other (2i?/7, 7: Lorentz factor) is ten times shorter at RHIC. Therefore the 
"pre-resonance" cc states collide with passing-by nucleons shortly after their pair 
production, when they are still small in size (not far from their production vertex) 
and thus have a smaller dissociation cross section compared to at SPS energy. 
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Fig. 36. J/tp Raa (normalized to Drell-Yan pairs) vs. L for various p+A and S+U 
systems. Note that this "L" is the effective average length travelled by (pre-) 
charmonium states, not by the initial state partons. This is different from the 
"L" in Fig. |35l The figure is taken from Ref . |105] . 

3. Implementation of CNM Effects 

In our transport approach charmonium production in the pre-equilibrium stage serves 
as the initial condition for the evolution in the hot medium. For simplicity we assume 
that the initial phase space distribution of charmonia, f{x,p,To), determined from 
\E' production in the pre-equilibrium stage, can be factorized into coordinate and 
momentum spaces. According to boost-invariance we only need to consider the phase 
space distribution in the transverse plane. 



/m> (6, ft, Pi, To) = U{b,Xt,To)U{b,Pt,TQ) 



(5.13) 



Let us start with the spatial part, fiii{b,Xt,To)- The Glauber model tells us the 
initially produced charmonia distribution according to the nuclear overlap function, 
f<ii{b,Xt,To) oc TA{xt)TB{xt — b). We augment it with the nuclear absorption: The 
probability for the dissociation of a pre-resonance cc state created at {xi.za) inside a 
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Fig. 37. Centrality dependence of J/ip RdAu in 200AGeV d+Au collisions measured 
by PHENIX [TnH] . 
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Fig. 38. Schematic representation of the geometry of J /if) nuclear absorption. The 
dot represents a (pre)-charmonium. The dashed arrows are their absorption 
trajectories. Figure taken from Ref. |99j . 

nucleus A through passing-by nucleons in the same nucleus A is 

oo 

Plsi.^t,zx) = (Ts.hsTA>{xuZA) with TA>{xt, za) = PA{xt, z) dz , (5.14) 

ZA 



where p is normalized to 1, see Eq. (15. 6p . The effective absorption cross section o"abs 
parameterizes the inelastic scattering between \1/ and nucleons. The integration limits 
are determined by the absorption trajectories, as illustrated in Fig. [381 Therefore the 
corresponding survival probability is 



vA-l 



P^^^^{Xt, Za) = {I - (7e,i,sTA>{Xt, Za)) ^ CXp (-(A - 1) (Tabs Ta> (f t, 2;a)) , (5.15) 

where the factor of (A—l) reflects the fact that the cc producing nucleon does not 
participate in the absorption. Similarly the probability for a pre-resonance cc state 
created at {xt — 6, zb) surviving from collisions with nucleons in nucleus B is 



-Psurv(^i - b, Zb) ^ exp -(5 - 1) (Tabs 7B>(ft - b, Zb] 



(5.16) 
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Therefore the number of surviving charmonia at Xt for impact parameter b is 

^'^pp f 

f^{b, Xt, To) =Ay—^AB / dzA dzBPA{xt, za) Psixt - b, zb) 

X PLvi^t, ZA)Psurv(^t - b, Zb) 

=Ay—^AB / dzA dzBpA{xt, za) psixt - b, zb) 

X exp <^ -(A - 1) / dzpA{xt, z)aabs ? 

X exp < -{B - 1) / dz'pB^Xt - 6, z!)aahs \ , (5-17) 

where a* is the charmonium production cross section in elementary nucleon-nucleon 
colhsions. The total number of surviving charmonia is thus 

N^{b)= f h{b,Xt,To)d^Xt. (5.18) 

It is convenient to define a nuclear suppression factor, 

Snuc{b) = TTT / d^XtdzAdzBpAiXuZA) pB{b-Xt,ZB) 

-lab{0) J 
X exp <^ -(A - 1) / dzpA{xt, z)aabs > 

X exp <^ -(5 - 1) / dz'psixt - b, z')aabs > , (5.19) 

which can be used to express the number of surviving charmonia as 

N^{b) = Ay^ ABTABib)S^ucib) . (5.20) 

dy 

Eqs. (l5.17p . (I5.18p . (I5.19P and (I5.20p can be reduced to p+A colhsions by setting 
PB{xt, z)=5^'^\xt)6{z) and B=l. For example, for p+A collisions, Eq. (I5.19P reduces 
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to 



Snuc{b) = — — / dzA pAib, Za) 

X exp <^ -{A - 1) / dzpA{b, z)aabs [ ■ (5.21) 

We parameterize both nuclear shadowing and nuclear absorption with the ef- 
fective \1/-N absorption cross section, fXabs- Applying Eq. fl5.2ip to p+A collisions at 
SPS we obtain a^l^ =7.3±lmb from the recent NA60 data at _Eiab=158GeV (corre- 
sponding to ^/sNN=^7.3GeV) |106] . This updated measurement at 158 GeV gives 
a significantly larger value than previously available for 400 GeV proton projectiles, 
aabs ^4.4mb [M] (the latter has been confirmed by NA60 [106], i.e., at 400 GeV), 
which has been used in our previous calculations [lU 1108] . The comparison with 
recent PHENIX data [I03l |I09] yields cXabs ^ 3.5 mb (5.5 mb) for v^=200AGeV Au- 
Au colhsions at mid rapidity, \y\ < 0.35 (forward rapidity, \y\ G [1.2,2.2]), see, e.g., 
Fig. [391 For simplicity, we assume the same absorption cross sections for the Xc as for 
the J/ip. However, for excited states dabs is expected to be significantly larger, even 
if they are not fully formed when the dissociation occurs. Taking guidance from the 
NA50 measurement with 400 GeV protons, we use a^^^ ~ 13 mb at y^=17.3AGeV 
and a^^^ ~ 6.5(10) mb at A/i=200 AGeV for mid (forward) y. 

The rather pronounced rapidity dependence of (Tabs at RHIC casts doubt on inter- 
preting this quantity as an actual absorption cross section. It seems more reasonable 
to associate its increase at forward y with nuclear shadowing [110] since the dissoci- 
ation kinematics is very similar between mid and forward rapidity. While this does 
not affect the use of our "effective" (Tabs; it does imply a nuclear shadowing effect on 
the open-charm cross section in A+A collisions (which is an important ingredient in 
the calculation of regeneration). As a "minimal" scheme we therefore associate the 
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Fig. 39. Centrality dependence of J/^ RdAu in 200AGeV d+Au collisions measured 
by PHENIX [103j together with a Glauber fit with a nuclear absorption cross 
section craf,s= 3. 2mb, 3.5mb, and 7.8mb for backward, mid- and forward rapid- 
ity. 

additional absorption of the J/ip yield at forward y (relative to mid rapidity) with 
a suppression of open charm production caused by shadowing, while we assume no 
shadowing corrections at mid rapidity. Thus, at both SPS and RHIC the number 
of primordially produced cc pairs at mid-rapidity is calculated from the p+p cross 
section as 

N^^'^ib) = Ay-^ ABTABib) , (5.22) 

while for forward y at RHIC we use 



Nl-{b) = Ay '^ 



dy 



y=1.7 



Cfor 



(5.23) 



Here, T^sib) is the usual nuclear overlap function, Eq. (15. 7p and S^^c, defined in 
Eq. (I5.19p . denotes the J/ip suppression factor due to CNM effects, parameterized by 
(Tabs in the Glauber formula, Eq. (I5.17p . In particular, the ratio S^^^/S^^ represents 
the extra suppression associated with nuclear shadowing, operative for both J/ip and 
cc production. 

For the momentum dependent part /*(6, Xt, tq) we take the charmonium pt spec- 
tra in p+p collisions as the baseline and apply a pt broadening associated with the 
Cronin effect. The J/ip pt spectra in p+p collisions at SPS are parameterized with 
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an exponential distribution, 

^ oc ^e-?/<^?) , (5.24) 

where (p^)=1.15GeV^ |104j . At full RHIC energy the J /ip pt spectra are parameter- 
ized with a power-law distribution, 

'' oc--^--^, (5.25) 



2TTptdpt [1 + {pt/oy]^ 

where D is adjusted so that (p^)=4.14 (3.56)GeV^ for J/^p produced at mid (forward) 
rapidity [95] and C is a normalization factor. 

The Cronin effect is readily implemented into the 3-momentum dependent part 
fmih.pt, tq) via a Gaussian smearing of the charmonium pt distribution in p+p colli- 
sions, ff{pt), 

The nuclear increase of the average pj^ {^Vt) — {vDaa — {vDpp^ i^ estimated within a 
random-walk treatment of parton-nucleon collisions [lllj as being proportional to the 
mean parton path length, (/"''), in the cold medium: (Apf) = agN {1°"^). The coefficient 
QgN is estimated from p+A data at SPS [104J and d-Au data at RHIC |103j . We use 
0^^=0.076 GeVVfm for v^=17.3AGeV Pb-Pb collisions and ag7v=0. 1(0.2) GeVVfm 
for 1/5=200 AGeV Au-Au collisions at mid (forward) rapidity. The mean parton path 
length (/°^) is determined from the geometric lengths /" and f which two partons a 
and h travel before they collide, weighted according to the survival probability of final 
state charmonia |lllj 

/ (fxt dZa dZh {l"-{Xu Za) + fih- Xu Zb)) K{b, Xt, Za, Zb) 

r\h) = ^ ^ ^ , (5.27) 

/ d^Xt dZa dZb K{b, fj, Za, Zb) 
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where 



t{Xt,Za) = A dzpAiXt,z)/po , 

^— oo 

l\b- Xt,Zb) = B / dzpsixt - b, z)/po 



(5.2J 



26 



are the respective path lengths travelled by parton a and b weighted according to the 
nuclear density p{s, z). The kernel K{b, Xt, za, zb) is given by 

K{b,Xt, Za, Za) = PA{Xt, Za)pB{Xt " &, Zb) 



X exp ( -aabs 

-26 



{A -I) / dzpA{xt,z) 



+ {B-1) dzpBixt-b,z) 

J —oo 

reflecting the survival probability of final state charmonia, recall Eq. (]5.17p . 



(5.29) 
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CHAPTER VI 

APPLICATION TO HEAVY-ION COLLISIONS 
In previous chapters we have introduced all the components of the kinetic rate equa- 
tion/transport approach. In this section we elaborate its application to heavy-ion 
collisions and compare our numerical results with experimental data. 

In Section IVIlRl we discuss the procedure of applying the kinetic approach to 
heavy-ion collisions. For the suppression of primordial charmonium {^=J /ip^ Xc-, 4'') 
we employ the Boltzmann equation to evaluate the time evolution of the \E' phase 
space distribution function. To solve the time evolution of regenerated "$ using the 
Boltzmann equation requires calculating \E' regeneration rates from the full time- 
dependent phase distribution function of charm quarks, which is numerically rather 
involved and is still a work in progress. For the most parts we adopt a simplified 
procedure: We estimate the inclusive yield from regeneration with a rate equation 
and estimate their p^-spectra with the blastwave formula. 

In Section IVIlB] we present the numerical results of the inclusive J/ip yield and 
its transverse momentum (pt) spectra and compare with SPS and RHIC data. We 
specifically compare two scenarios with the internal ( "strong binding" ) or free ( "weak 
binding") energy identified as the heavy quark 2-body potential. The effects specifi- 
cally relevant for high pt (>5GeV) J/ip production are discussed. The discussions in 
Section EHX] and Section ElB] mostly follow Ref. [2B]. 

In Section IVIlCl we study the impact on charmonium regeneration due to off- 
equilibrium effects in charm-quark phase space distributions. Specifically, we com- 
pare the inclusive yield and pt spectra of regenerated J/ip from limiting charm-quark 
spectra including: (1) thermal charm-quark spectra; (2) pQCD charm spectra. Also, 
the impact of a charm-quark correlation volume is discussed within the framework of 



98 



the Boltzmann equation. 

In Section |VI1[D] the production of excited \l/ states such as Xc and ip' at SPS and 
RHIC is briefly discussed. 

Finally, in Section IVIIlEJ we present predictions for charmonium production at 
FAIR energies, where the medium is expected to have a lower initial temperature and 
higher baryon density compared to SPS energies. 

A. The Rate Equation Approach 

Throughout this chapter we solve kinetic equations (Boltzmann or rate equation) 
separately for '^=J /ip, Xci '4^' ■ In the hot medium the dissociation and regeneration 
are the two main processes affecting the \E' yield. It is often desirable to disentangle 
these two effects and study their respective strength. For this purpose we decompose 
the charmonium distribution in the medium at any time r, 

/vt(pt, xt, r) = fl"'^{pu xu r) + f^^^{pu Xt, r) , (6.1) 

into a (suppressed) primordial component and a regenerated one by exploiting the 
linearity of the Boltzmann or rate equation. According to boost-invariance we only 
need to solve for the time evolution of f^{pt, Xt, t) in the transverse plane. We define 
f^^^^iPt, Xt, t) as the solution of the homogeneous Boltzmann equation, 

dft-^/dr + vt ■ Vjf^ = -a^ff"^ , (6.2) 

with the same initial condition as for the full Boltzmann equation, (12. ip . f^^^ipt, Xt, tq) 
= fqi{pt,Xt,TQ), which is obtained from Eq. (I5.13p . The explicit expression for the 
solution of Eq. (16. 2p is 

ff'^iPt, Xt, t) = Mpt, Xt -Vt{r- To), ro)e' A '^^'-*b*----(-^'),^') _ (g^g) 
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We account for the "leakage effect", i.e., charmonia escaping tlie fireball volume no 
longer being subject to suppression, by setting a^ = whenever \xt — Vtir — t')\ > 
R{t'), where -R(t') is the fireball radius at time r'. Due to the leakage effect the 
suppression of high pt charmonia is reduced compared to low pt ones, since the former 
are more likely to escape from the fireball. 

The regeneration component, f^^, follows as the difference between the solution 
of the full and the homogeneous Boltzmann equation, which can be expressed as 

d/r/dr + vt ■ Vjr = -a'^ir + P^ , (6.4) 

with vanishing initial condition, f^^{x,p, t < r*) = 0. The onset time of regeneration 
processes, r*, is defined by the dissociation temperature T{tq) = T^^^^ for each state 
\1/. The explicit expression for the solution of Eq. (16 ■4p is 

Due to the complication mentioned in the introduction of this chapter we adopt 
the following approximation: Instead of the Boltzmann transport equation (16.41) we 
solve the following rate equation for the inclusive yield of the regeneration component, 

^ = -r^ {Nl^^ - N'^) . (6.6) 

For the dissociation rate, F^, we employ a 3-momentum average whose precise value is 
obtained by matching the final yield of the loss term to the exact result obtained from 
solving the momentum- dependent Boltzmann equation for the primordial component, 
Eq. (16. 2p . For the \1/ equilibrium limit we adopt the statistical model mentioned in 
Section HTIICl Since the thermal production and annihilation rates of cc are believed to 
be small at SPS and RHIC energies, cc pairs are assumed to be exclusively produced 
in primordial N+N collisions and conserved thereafter. The open and hidden charm 
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states are then populated in relative chemical equilibrium according to the canonical 
charm-conservation equation, (12 .49 p . 

iVe, = Zi^i|||i| + iVhid, (6.7) 

with Ncc'- total number of charm-quark pairs from initial production; Zi = 7cVfb^op: 
1-body open charm partition function (recall Eq. f l2.38p ). with pertinent equilibrium 
density nop; N^id = Tc^B'^hid: total number of all charmonium states with pertinent 
equilibrium density rihid; and jc'- charm-quark fugacity accounting for the deviation 
of chemical equilibrium with the heatbath (7c=l in full equilibrium). The ratio of 
modified Bessel functions, Ii{2Zi)/Iq{2Zi), on the right-hand-side of Eq. (16. 7p is the 
characteristic canonical suppression factor which accounts for the exact conserva- 
tion of net-charm number, Nc — Nc, in each event |112| 147] : for Zi <^ 1, one has 
Ii{2Zi)/Iq{2Zi) —7- Zi, which acts as an additional (small) probability to enforce a 
vanishing net charm content in the system [i.e., both c and c have to be present 
simultaneously), see Section HTIIClIT] for a more detailed discussion. 

The 1-body open charm partition function, Zi, is evaluated as follows. For the 
QGP phase in the weak-binding scenario only charm quarks are counted as open- 
charm states. In the strong-binding scenario, the T-matrix calculations of Ref. [50] 
suggest that cq and cq (charm-light) bound states (D-mesons) survive in QGP up to 
~1.3 Tc] therefore, we count both charm quarks and the lowest-lying 5- wave D-mesons 
{D, D*, Dg and D*), as open charm states for T < 1.3Tc. The charm-quark masses in 
the QGP correspond to the temperature-dependent ones displayed in Fig. [161 while 
for the meson resonances above T^ we estimate from Ref. [50] rriu = m£,* ~ 2.0 GeV 
and m£,^ = m£,* ~ 2.1 GeV (hyperfine splitting has been neglected). For the hadronic 
phase all charmed hadrons listed by the particle data group [113] are counted as open- 
charm states, with their vacuum masses. The number of hidden charm states, A^hid, 
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Fig. 40. Temperature dependence of the charm-quark fugacity (left panel) and the 
in-medium J /ip equilibrium limit (right panel) using the statistical model 
in the QGP within the strong-binding scenario (dot-dashed lines: with and 
without D-meson resonances below and above 1.3Tc~234MeV, respectively; 
solid line: smooth interpolation of the previous two cases; see text for de- 
tails), the weak-binding scenario (dashed line) and in the hadronic phase for 
temperatures below Tc=180MeV. 

is evaluated in line with the existing charmonium states and their masses at given 
temperature T, but its contribution to N^c is numerically negligible. 

Knowing riop, "n-hid and Vyb at each temperature, one can solve Eq. (16.71) for the 
charm-quark fugacity, 7c(7'), and apply it to compute the statistical equilibrium limit 
of each charmonium state as 



m = 7c ^FB n^ , 



(6.8) 



in terms of its equilibrium density, riij,. In Fig. |10]we collect the numerical results of 
the charm-quark fugacity 'ydT) and the statistical equilibrium limit for J/ip abun- 
dances (excluding feeddown) for central 200AGeV Au+Au collisions at RHIC. The 
discontinuity at 1.3 Tc for the strong-binding scenario (dot-dashed line) is due to 
the inclusion of the D resonances in the QGP medium. We smoothly interpolate 
around the melting temperature for the D-mesons with a hyperbolic tangent func- 



102 



tion (solid line) to represent a more gradual (dis) appearance of the D resonances (we 
have checked that this procedure has negligible impact on the calculation of observ- 
ablesin Section IVHH]). 

To achieve a more realistic implementation of the statistical equilibrium limit, we 
apply two corrections to N^'^'^ to schematically implement off-equilibrium effects of 
charm quarks in momentum and coordinate space. The former is aimed at simulating 
incomplete thermalization of the charm-quark pt spectra throughout the course of the 
thermally evolving bulk medium. We have shown in Section HTTllCl that the coalescence 
rate from non- or partially thermalized c- and c-quark spectra is smaller than for fully 
thermalized ones [114[ [58] , since the former are harder than the latter and thus provide 
less phase-space overlap for charmonium bound-state formation. We implement this 
correction by multiplying the charmonium abundances from the statistical model with 
a schematic relaxation factor [114J, 

N'^ = n{T) Nl"'' , 7^(r) = 1 - exp(-r/r:'^) , (6.9) 

where r^^ is a parameter which qualitatively represents the thermal relaxation time of 
charm quarks (it is one of our 2 main adjustable parameters in our phenomenological 
applications in Section IVllBj) . A rough estimate of this time scale may be obtained 



from microscopic calculations of this quantity within the same T-matrix approach 
as used here for charmonia, where the thermal charm-quark relaxation time turns 
out to be Tgq ~3-10fm/c |115[ [50| . Such values allow for a fair description of open 
heavy-flavor suppression and elliptic flow at RHIC |115[ 1116] . The second correction 
is applied in coordinate space, based on the realization that, after their pointlike 
production in hard N-|-N collisions, the c and c quarks only have a limited time to 
diffuse throughout the fireball volume. At RHIC and especially at SPS only few cc 
pairs are produced {e.g., dNcs/dy ~ 1.2 in semicentral (6=7 fm) Au+Au collisions 



103 



at RHIC), and the hadronization time (QGP lifetime) is smaller than the fireball 
radius. Thus, c and c will not be able to explore the full fireball volume but rather 
be restricted to a "correlation volume", \4o [S] (the analogous concept has been 
successfully applied to strangeness production in p+A and A+A collisions in the SPS 
energy regime [117j ). see Section [IT1IC1I2] for a more detailed discussion. In the present 
rate equation approach we implement this correction by replacing the fireball volume 
VpB in the argument of the Bessel functions in Eq. (16. 7p by the correlation volume 
Ko l^lllSj . The latter is identified with the volume spanned by a receding cc pair, 

4:71 

Ko(r) = y(ro + (t;e)r)3, (6.10) 

where tq ~ 1.2 fm represents an initial radius characterizing the range of strong 
interactions, and (vc) is an average speed with which the produced c and c quark 
recede from the production point; we estimate it from the average pt in D-meson 
spectra in p+A colhsions [Il9l Hi |55] as (vc) ~ 0.55(0.6)c at SPS (RHIC). The 
correlation volume leads to a significant increase of •jc (since lo/h is reduced) and 
thus of the modified \Ef "equilibrium limit" due to locally increased cc densities. 

The Pt spectra of regenerated charmonia are approximated by local thermal dis- 
tributions boosted by the transverse fiow of the medium, amounting to a standard 
blastwave description [79], Eq. (I4.14p . 

-^«m,y^ ,*/<■, (^^j /„ (^j (6,11) 



(rrit = \/vr?^ + pW The medium is characterized by the transverse- fiow rapidity 
p = tanh^^ vjf\r) using a linear fiow profile ft(r) = Vg^ with a surface velocity v^ = 
oxTmix and transverse fireball radius R = R{T^ix) as given by the fireball expansion 
formula, Eq. (14. ip . at the end of the mixed phase, Tmix- We evaluate the blastwave 
expression at the hadronization transition [Tc) and neglect rescattering of \I''s in the 
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hadronic phase. Two additional effects are neglected in this treatment, which, to a 
certain extent, tend to compensate each other: on the one hand, due to incomplete 
charm-quark thermalization, one expects the regenerated charmonium spectra to be 
harder than in equilibrium, but, on the other hand, a good part of the regeneration 
occurs before the mixed phase |120l BO] (as will be shown in Section IVIlBJIT]) so that 
the evaluation of the blastwave expression at the end of the mixed phase presumably 
overestimates the blue shift due to the flow field. We will explicitly check these effects 
in Section IVIlCl Ultimately, an explicit evaluation of the gain term with realistic 
(time-dependent) charm-quark spectra within a Boltzmann equation |121] will be 
able to lift these approximations. 

B. Inclusive J/^ Yield and pt Spectra at SPS and RHIC 

In this section we present and discuss the numerical applications of the above frame- 
work to J /ip data in URHICs at SPS and RHIC. For each observable, we confront the 
results of the strong- and weak-binding scenario in an attempt to discriminate qual- 
itative features. The feeddown to J /il) from Xc and ip' states is taken into account, 
assuming fractions of 32% and 8%, respectively, for primordial production in pp colli- 
sions. We have divided the discussion into the centrality dependence of inclusive J /ip 
yields in Section IVlJiBJIT] and the pt dependence of J/iIj in Section IVlJlBJI^l 

1. Inclusive J /ip Yield 

The J /tp yield in A+A collisions is usually quantified in terms of the nuclear modifi- 
cation factor as a function of centrality. 
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Fig. 41. J/ip production (normalized to Drell-Yan pairs) versus centrality at SPS 
evaluated with the thermal rate-equation approach, compared to NA50 
data |122l 1123] . Solid lines: total J/ip yield; dashed lines: suppressed pri- 
mordial production; dot-dashed hues: regeneration component; dotted lines: 
primordial production with CNM effects only. Left panel: strong-binding 
scenario; right panel: weak-binding scenario. 

where A^coii(^) is the number of binary collisions of the incoming nucleons at impact 
parameter b, recall Eq. (15. Sp . Before we turn to the results, we recall the two main 
parameters in our approach, which are the strong coupling constant, as, and the 
thermal charm-quark relaxation time, r^^. The former controls the inelastic charmo- 
nium reaction rate, (Eqs. (I3.10p and (I3.15P ). and the latter the magnitude of the \Ef 
equilibrium limits, Eq. (I6.9p . We adjust them to reproduce the inclusive J/ip yield 
for central A+A collisions at SPS and RHIC, within reasonable bounds. For a^ we 
find that a common value of 0.32, which is at the upper end of the value in the 
Coulomb term in the QQ free energy, can be used, in combination with r^'^=3.8fm/c 
for the strong-binding scenario and r^'^=1.6fm/c for the weak-binding scenario. For 
simplicity, we refrain from introducing an additional temperature dependence into 
these parameters. The composition of the total yield, its centrality dependence and 
Pt spectra can then be considered as a prediction within each of the 2 scenarios. 
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Fig. 42. J ftp Raa versus centrality at mid-rapidity at RHIC evaluated with the ther- 
mal rate-equation approach, compared to PHENIX data [124J. Solid lines: 
total J/ip yield; dashed lines: suppressed primordial production; dot-dashed 
lines: regeneration component; dotted lines: primordial production with CNM 
effects only. Left panel: strong-binding scenario; right panel: weak-binding 
scenario. 

We begin with J/ip production in yi=17.3AGeV Pb+Pb collisions at SPS, for 
which we compare our results in the strong- and weak-binding scenario with NA50 
data in Fig. |4T1 For these data, the denominator in Eq. (I6.12P is replaced by the num- 
ber of Drell-Yan dileptons at high mass, while the numerator includes the branching 
ratio into dimuons. The pertinent proportionality factor, equivalent to the pp limit of 
this ratio (47. Oil. 4 |125j ). and the CNM-induced suppression (dotted line in Fig. HI]) 
are inferred from the latest NA60 p+A measurements |106j . which we reproduce us- 
ing the Glauber model formula, Eq. (I5.17p . with (Tabs=7.3mb. The suppression of 
the primordial component (dashed line) relative to nuclear absorption (dotted line) 
represents the "anomalous" suppression by the hot medium, which increases with 
centrality due to higher initial temperatures and longer fireball lifetimes. The regen- 
eration component increases with centrality as well, mostly due to the increase of 
the 7^-factor and the larger lifetime which facilitates the approach to the equilibrium 
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Fig. 43. J/i/j Raa vs. centrality at forward rapidity compared to PHENIX data |124] . 
Solid lines: total J/ip yield; dashed lines: suppressed primordial production; 
dot-dashed lines: regeneration component; dotted lines: primordial produc- 
tion with CNM effects only. Left panel: strong-binding scenario; right panel: 
weak-binding scenario. 



limit according to Eq. f l6.6p . Because of detailed balance between dissociation and 
regeneration, an increase in the former also implies an increase in the latter. The 
sum (solid line) of primordial and regeneration contributions describes the centrality 
dependence of the inclusive J/ip yield at SPS reasonably well in both scenarios. In 
the strong-binding scenario the primordial component is dominant and the majority 
of the anomalous suppression originates from the dissociation of Xc and ip', since at 
the temperatures realized at SPS (Tq ~ 200 MeV) the quasifree dissociation rates for 
Xc and ip' are much larger than those for J/ip, recall Fig. |2T1 In the weak-binding sce- 
nario, however, the regeneration yield becomes comparable to the primordial one for 
semi-/central collisions due to larger dissociation rates and the smaller charm-quark 
equilibration time scale. 

Next we examine the centrality dependence of J/ip production in 200AGeV 
Au+Au collisions at RHIC, first focusing on mid rapidity {\y\ < 0.35), as shown 
in Fig. US The suppression due to CNM effects (dotted line in Fig. H2l) is inferred 
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Fig. 44. Ratio of Raa for J/ip at forward and mid rapidity versus centrality in the 
strong (solid line) and weak (dashed line) binding scenarios compared to 
PHENIX data |124j . In the upper two curves, CNM effects have been di- 
vided out in both numerator and denominator of the ratio. 

from latest PHENIX d+Au measurements, which we reproduce using the Glauber 
model formula Eq. (I5.17P with o"abs=3.5mb. For A^part — — 100, the composition of 
primordial and regeneration contributions is quite comparable to the SPS for A^part — 
— 400 within both scenarios. Beyond Apart — 100, suppression and regeneration 
continue to increase, leveling off at an approximately 50-50% (20-80%) partition for 
primordial and regeneration in the strong-binding (weak-binding) scenario in central 
collisions. 

Let us now turn to J/ip production at forward rapidity {\y\ G [1.2, 2.2]) at RHIC, 
shown in Fig. HHl Again, both strong- and weak-binding scenarios reproduce the 
experimental data fairly well, with similar relative partitions for primordial and re- 
generation contributions as at mid rapidity. However, one of the "puzzles" about 
J/ip production at RHIC is the fact that the total J/ip yield is more strongly sup- 
pressed at forward rapidity than at mid rapidity. In our approach this follows from 
the stronger shadowing at forward rapidity leading to less primordial production for 
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both J /%l) and cc pairs. The former (latter) leads to a reduction of the primordial 
(regeneration) component. Since the thermodynamic properties of the fireball are 
quite similar at mid and forward rapidity (recall Fig. I29l) . charmonium suppression 
and regeneration in the hot medium are very similar between the mid and forward 
rapidity as discussed in Ref. [76] . To quantify the difference at forward rapidity and 
mid rapidity we display the ratio between the corresponding Raa^ in Fig- SH for 
both scenarios, which clearly illustrates the importance of CNM effects to properly 
reproduce the data. 

In Section IIIIIlA] we have argued that the strong- and weak-binding scenarios 
discussed here may be considered as limiting cases for J /ip binding in the QGP, as 
bracketed by the identification of the heavy-quark internal and free energies with a 
Q-Q potential. From the results above we believe that these scenarios also provide 
a reasonably model-independent bracket on the role of suppression and regeneration 
effects, in the following sense: At SPS, the strong-binding scenario defines a "mini- 
mal" amount of dissociation required to provide the anomalous suppression beyond 
CNM effects (a small regeneration component is inevitable due to detailed balance). 
The application to RHIC energy then implies an approximately equal partition of 
primordial and regenerated charmonia in central Au+Au, not unlike Ref. ^40j where 
the vacuum charmonium binding energies ("strong binding") have been used in the 
QGP (together with the gluo-dissociation). On the other hand, in the weak-binding 
scenario, a large part of the J /ip yield in central A+A is due to regeneration even 
at SPS, limited by the constraint that for sufficiently peripheral collisions (and at 
sufficiently large pt) a transition to primordial production compatible with p+A data 
should be restored. Clearly, for central A+A at RHIC (and certainly at LHC) the final 
yield is then dominated by regeneration. Since both scenarios describe the inclusive 
yields reasonably well, it is mandatory to investigate more differential observables to 
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Fig. 45. J/ip abundance as a function of time in central (A^part=380) Au+Au colli- 
sions at RHIC. Solid lines: total; dashed line: primordial component; dot- 
dashed line: regeneration component; double dot-dashed line: equilibrium 
limit, NjJ,. Left panel: strong-binding scenario; right panel: weak-binding 
scenario. 

find discriminating evidence. This will be pursued in the following section. 

It is instructive to examine the time evolution of J/ip production in the two 
scenarios, displayed in Fig.|l5]for central collisions at mid rapidity at RHIC (excluding 
feeddown from Xc and ip'). In both scenarios most of the dissociation and regeneration 
indeed occur in the QGP and mixed phase, since the hadronic reaction rates are small. 
In the weak-binding scenario the time-dependent J/ip yield exhibits a "dip" structure 
around r ^ 1.5 fm/c because the large dissociation rates suppress primordial J/ip very 
rapidly and regeneration only starts after the medium temperature falls below the 
J/ip dissociation temperature {Tjff ~ 1.25 T^). This scenario is closest in spirit to 
the statistical hadronization model [36] where all initial charmonia are suppressed (or 
never form to begin with, except for corona effects) and are then produced at the 
hadronization transition. 
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2. J /ip Transverse Momentum Spectra 

a. Average Transverse Momentum 

The results of the previous section suggest that, within the current theoretical {e.g., 
charm-quark relaxation time, r^"^) and experimental uncertainties both of the "lim- 
iting" scenarios can reproduce the centrality dependence of the inclusive i?yj^(A^part) 
reasonably well at both SPS and RHIC energies. However, the composition between 
suppression and regeneration yields is rather different which ought to provide a key 
to distinguish the two scenarios. The obvious "lever arm" are charmonium pt spec- 
tra [37]. One expects that the primordial component is characterized by harder 
Pt spectra (following a power law at high pt) while the regeneration component pro- 
duces softer Pt spectra characterized by phase-space overlap of (partially) thermalized 
charm-quark spectra. However, in practice, the transition from the "soft" recombina- 
tion regime to the "hard" primordial regime is quite uncertain; e.g., collective flow and 
incomplete thermalization of c-quarks can lead to a significant hardening of the regen- 
erated J /ip spectra, while a dissociation rate which increases with 3-momentum |108j 
can induce a softening of the spectra of surviving primordial charmonia. 

For a more concise discussion of the pt dependence of J /ip as a function of cen- 
trality at SPS and RHIC we here focus on the average pi, as compiled in Figs. HE] 
andHTl At the SPS (Fig. |16]), the centrality dependence of {pi) is largely dictated by 
the the Cronin effect in the primordial component, especially in the strong-binding 
scenario where this contribution dominates the yield at all centralities. The momen- 
tum dependence of the dissociation rate induces a slight suppression of {pD at large 
centrality compared to the case where only CNM effects are included (dashed vs. 
dotted line) |108j . In the weak-binding scenario, larger contributions from regenera- 
tion induce a slight "dip" structure at intermediate centralities due to a rather small 
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Fig. 46. {p"^) of J/ip vs. centrality at SPS, compared to NA50 data [1261 [IDl]. In each 
panel, (pf) is plotted for total J/ip yield (primordial + regeneration com- 
ponent; solid lines), the suppressed primordial component (dashed line), the 
regeneration component (dash-dotted line) and primordial production with 
CNM effects only (dotted lines). The left (right) panels correspond to the 
strong-binding (weak-binding) scenario. The transverse energy, Et, is a mea- 
sure of centrality. Its value in GeV is about 0.275 xA'paj-t- 

collective flow at the end of the (relatively short) mixed phase in these collisions. 
At RHIC energy (Fig. W7}i . the individual primordial and regeneration components 
show qualitatively similar behavior for (Pi)(A^part) as at SPS, i.e., an increase due to 
Cronin effect and collective flow, respectively. At mid rapidity, the general trend is 
that with increasing centrality the growing regeneration contribution pulls down the 
average {pf), in qualitative agreement with the data. The curvature of the (p^)(A^part) 
dependence, which appears to be negative in the data, is not well reproduced, nei- 
ther by the strong- nor by the weak-binding scenario, even though the deviations are 
smaller in the former. A more microscopic calculation of the gain term, together with 
more accurate estimates of the Cronin effect, are warranted to enable more definite 
conclusions. For both rapidity regions, the (p'f) of the suppressed primordial compo- 
nent is slightly larger in the weak- than in the strong-binding scenario. This is caused 
by the stronger 3-momentum increase of the dissociation rate in the strong-binding 
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Fig. 47. (p^) of J li\) vs. centrality at RHIC at mid-rapidity (upper panels) and forward 
rapidity (lower panels), compared to PHENIX data [124] . In each panel, (p^) 
is plotted for total J/?/' yield (primordial + regeneration component; solid 
lines), the suppressed primordial component (dashed line), the regeneration 
component (dash-dotted line) and primordial production with CNM effects 
only (dotted lines). The left (right) panels correspond to the strong-binding 
(weak-binding) scenario. 
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scenario, recall the right panel of Fig. [2TJ 

The overall comparison to SPS and RHIC data for the pt dependence of J/tp^s 
seems to indicate a slight preference for the strong-binding scenario. This is mostly 
derived from the observation that for peripheral collisions the experimentally ob- 
served (j9j) essentially follows the extrapolation of the Cronin effect, suggesting J/ip 
production of predominantly primordial origin (the collective flow imparted on the 
regeneration component appears to be too small at these centralities, rendering the 
weak-binding scenario problematic). 

b. High Pt J/ip Production and Elliptic Flow 

The average pf essentially characterizes the momentum dependence of charmonium 
production at low and moderate pt where most of the yield is concentrated. Recent 
RHIC data [92| |93] have triggered considerable interest in J ftp production at highp^ ~ 
5 — lOGeV which is expected to provide complementary information. It was found 
that the suppression in Rj^{pt > 5 GeV) in Cu-Cu collisions is reduced compared 
to the \ow-pt region, with i?^^-values of ~0.7-l or even larger [93]. This is quite 
surprising in light of the light-hadron spectra measured thus far at RHIC which all 
exhibit stronger suppression of Raa — 0.25 for Pf > 6 GeV (even electron spectra from 
open heavy flavor, i.e., charm and bottom decays). It also appears to be at variance 
with the thermal J/tp dissociation rates which, if anything, increase with momentum 
(recall Fig. [21]) and thus imply a stronger suppression at higher pf. Furthermore, 
the leakage effect mentioned in Section IVlllA] is not strong enough to produce the 
experimentally observed increase in RJ^{pt > 5 GeV) [108J. 

Therefore we consider the following two effects primarily relevant at high pt |127] : 
(1) Finite formation times for the charmonium states |128l 11291 I130J and (2) Bot- 
tom feeddown. Concerning (1), one expects reduced geometrical sizes for a "pre- 
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Fig. 48. Ratio of J/ip from 5-nieson feeddown to inclusive J/ip indicated by Teva- 
tron [m] and STAR data [T82] . 

resonance" cc pair relative to a fully formed charmonium due to a finite formation 
time, Tf, required to build up the hadronic wave function. If the "pre-resonance" cc 
pair is in a color-singlet state, as suggested by the color singlet model (CSM) (recall 
Section rvTRj) . its smaller geometrical size will translate into a smaller dissociation 
cross section, since it appears color-neutral for incoming partons whose momenta are 
not large enough to resolve its inner structure. For a schematic estimate we parame- 
terize the evolution of the pre-hadronic dissociation rate as 



rpre-*(T) = r^r/r' 



lab 
f 



Jab 
'/ 



r < rf = TfTTit/m^ 



(6.13) 



with r^: (nuclear, partonic or hadronic) dissociation rates for a formed charmonium, 
r: fireball proper time, rj=0. 89(2. 01, 1.50) fm/c: formation time of J/tp{xc, "^0 in its 
rest frame, and rrit = (m|, +PtY^'^. Essentially, (pre-) charmonium dissociation rates 
acquire an additional momentum dependence through Lorentz time dilation, being 
reduced at high pt. Also note that the longer formation times of Xc and ip' imply less 
suppression relative to J/ip, quite contrary to standard dissociation and regeneration 
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Fig. 49. J/ip Raa versus pt in central 200AGeV Au+Au collisions including forma- 
tion-time effects (fte) and i?-meson feeddown (Bfd) contributions. PHENIX 
data |124] are compared to our rate-equation calculations in the strong- and 
weak-binding scenarios (left and right panels, respectively). Upper panels: 
mid-rapidity; lower panels: forward rapidity. 



mechanisms: as higher excited cc states than J /if) they have smaller binding energy 
and therefore are more easily destroyed, and they are heavier so that their equilibrium 
abundances are suppressed compared to J /ip by the thermal Boltzmann factor leading 
to a reduced regeneration. Therefore high-p^ Xc and ip' could provide a rather unique 
signature of the formation time effect. Concerning (2), Fig. HH] shows recent data 
on the B — t- J /ip feeddown fraction in elementary p+p(p) collisions, which is quite 
significant. As an estimate of this contribution, we use the Tevatron data [131] and 
replace the corresponding fraction of primordial component. 
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Fig. 50. J/tp Raa vs. pt in 0-60% Cu+Cu collisions, compared RHIC data [921 [T33]. 
The formation time effects and 5-meson feeddown (green band) are included. 
The partition between the primordial and regeneration components is based 
on the thermal rate-equation approach in Ref. |108j . 



Combining formation time effects and 5-meson feeddown we obtain results for 
Rjj^{pt), as displayed in Fig. HO] up to pt = 10 GeV. We find that the suppression is 
reduced to about 0.5 at the highest pt, compared to about 0.4 at low pt- This is similar 
to the moderate enhancement we found in the Cu-Cu case |127j . as shown in Fig. [501 
Surprisingly, the high-p^ suppression turns out to be very similar in both strong- 
and weak-binding scenarios, despite the fact that the high-pt yield is exclusively due 
to the primordial component whose strength is very different in the 2 scenarios at 
low Pt. The reason is the 3- momentum dependence of the dissociation rates, which 
become quite similar in the 2 scenarios at large 3-momentum: at p ~ 10 GeV, the 
difference in the energy-threshold due to binding energies of several 100 MeV becomes 
less relevant so that a collision with almost any thermal parton is energetic enough 
for dissociating the bound state. For J/ip at forward rapidity (shown in lower panels 
of Fig. HHj), we additionally include the effect that the nuclear shadowing (responsible 
for the extra CNM-induced suppression relative to mid-rapidity) decreases with the 
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Pt of the primordially produced J/ip. We assume that it ceases to exist at around 
pj=10GeV, recall the discussion in Section rVllB]|2]|al 

Let us also estimate the elliptic flow, f2(pt), of the J/ip, (recall Eq. fll.l2p ). which 
is hoped to be another good discriminator of primordial and regenerated production. 
For the former, a nonzero V2 is basically due to the path-length difference when 
traversing the azimuthally asymmetric fireball, typically not exceeding 2-3% |K Hll34] . 
For the latter, much larger values can be obtained if the coalescing charm quarks are 
close to thermalized |135[ 158] . However, as pointed out in Ref . |127] , the well-known 
mass effect suppresses the V2{pt) for heavy particles at pt < m; it is precisely in this 
momentum regime where the regeneration component is prominent. In Ref. |127] we 
estimate the total J/ip V2{pt) by combining the V2 from the primordial production 
computed in Refs. |10l 1134] and that from the regeneration computed in Refs. |135[ 
[58] with relative fractions determined from our thermal rate-equation approach. In 
Fig. 151] we show the resulting total J /ip V2{pt) for 20-40% central Au+Au collisions, 
where neither the formation time effect nor 5-meson feeddown is included. These 
two effects would further increase the fraction of the primordial components in the 
intermediate to high pt region leading to an even smaller total V2 in this region. Thus, 
we predict that in both strong- and weak-binding scenarios the total J /ip V2{pt) does 
not exceed ~5% at any pj, similar to what is found in Refs. |10l 1136] . The only 
alternative option we can envision are strong elastic interactions of the J /ip which 
are only conceivable in the strong binding scenario to avoid break-up in scattering off 
thermal partons |55j . 
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Fig. 51. J /ip V2{pt) for 20-40% central Au+Au collisions at RHIC Figure taken from 
Ref. |127] . The band represents the sum of uncertainties based on differ- 
ences between two independent "input" V2{pt) calculations for each of the two 
components. 

C Explicit Calculations of the Regeneration Component 

In this section we explicitly calculate J /if) regeneration from input charm-quark spec- 
tra using the Boltzmann equation (16 ■4p instead of the rate equation (16. 6 p to check 
the sensitivity of ^ regeneration to off-equilibrium effects in charm-quark spectra. 



1. Sensitivity of J /ip Regeneration to Charm-Quark pt Spectra 

We first compare the inclusive yield of J/ip^s regenerated from the three types of 
charm-quark momentum spectra introduced in Section UTTllCl 1) thermal charm-quark 
spectra; 2) 3-dimensionally isotropic pQCD spectra; 3) transversely pQCD + lon- 
gitudinally thermal spectra. Our calculation is performed at RHIC energy where 
the regeneration component takes up a significant fraction. For comparison pur- 
pose we only consider J/ip regeneration in QGP phase with the quasifree process 
z + c + c— j-iH- J/ip {i = g, q, q) within the strong binding scenario. Neither the 
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Fig. 52. J/?/' regeneration from different charm-quark spectra. Double-dot-dashed line: 
thermal charm-quark spectra. Dashed line: 3-dimensionally isotropic pQCD 
charm-quark spectra. Dot-dashed line: pQCD charm-quark spectra in trans- 
verse plane -|- thermal spectra in longitudinal direction. Solid line: pQCD 
charm-quark spectra in transverse plane with angular correlation between c 
and c + thermal spectra in longitudinal direction. 

canonical ensemble effect nor the correlation volume effect (see Section IIT][C|) are ap- 
plied; these two effects will be discussed in the next section. The resulting centrality 
dependence of inclusive yields of regenerated J/ip is compared in Fig. |52j Clearly, 
the thermal charm-quark spectra are most efficient in regenerating J/ip due to the 
large phase space overlap between c, c quarks and light partons in the medium. The 
regeneration from longitudinally thermal and transversely pQCD spectra is reduced 
to about half relative to the thermal spectra. The regeneration from 3-dimensionally 
isotropic pQCD charm-quark spectra is even reduced by a factor of 4. This result is 
similar to what is found in a quark coalescence model [58] where the number of J/ip 
coalesced from pQCD charm-quark spectra is smaller than that from thermal spectra 
by a factor of 3. However in Ref. [ID] a similar Boltzmann transport approach was 
employed for J/\E' regeneration with the (inverse) gluo- dissociation process as the re- 
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Fig. 53. (p^) of J/ip^s regenerated from different charm-quark spectra. Double-dot- 
dashed hne: thermal charm-quark spectra; solid line: 3-dimensional isotropic 
pQCD charm-quark spectra; dot-dashed line: longitudinal thermal + trans- 
verse pQCD charm-quark spectra; dotted line: longitudinal thermal + trans- 
verse pQCD charm-quark spectra with angular correlation between c and c; 
double dash-dotted line: the (pf) from the rate equation + blastwave treat- 
ment; dashed line: the (pf) of the primordial component. 

generation mechanism and it is found that the inclusive yield of regenerated J/ip from 
pQCD charm spectra is quite comparable with that from thermal spectra (within a 
~30% difference). Further investigations are needed to clarify the discrepancy. We 
also find that if the angular correlation between the transverse momentum of c and 
c is included, to reflect the back-to-back charm production in initial hard collisions, 
the regeneration is further reduced to about 1/10 of the regeneration from fully ther- 
malized charm-quark spectra, as a result of further reduction of the cc overlap in 
momentum space. 

Next we proceed to the (pf) of regenerated J/ip, displayed in Fig. [531 As ex- 
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pected, the {pf) of J /if) regenerated from pQCD charm-quark spectra is larger than 
that from thermal spectra. We note that the {pi) from thermal charm-quark spectra 
is lower than that estimated from the blastwave formula. This is due to the fact that 
in the strong-binding scenario most of the regeneration processes occur at an early 
stage of the medium evolution (see Fig. HSj) . when the collective flow has not yet fully 
built up. Here the (transverse) medium flow effect is included only for the thermal 
charm-quark spectra, for the other two types of pQCD spectra we assume the medium 
is at rest in the lab frame. We also note that the introduction of the angular correla- 
tion between c and c momenta significantly lowers the {pf) of regenerated J /ip. This 
is because high pt J/ip^s are more likely to be regenerated from two charm quarks 
with a small angle between them, the probability of which is significantly suppressed 
by the back-to-back correlation of charm pairs. Therefore the regeneration of J/^/^'s 
at high Pt is more reduced than at low pt leading to a reduction of (pf). 

2. Sensitivity of J/tp Regeneration to Charm-Quark Correlation Volume 

In this section we explicitly evaluate J/ip regeneration from individual events with 
integer numbers of cc pairs produced, recall Fig. [T3 The correlation volume effect 
(see Section [ITEE]) is applied separately to events with integer (/c) charm quark pairs. 
At the end we average the regenerated J/tp over these events. The time dependence 
of the correlation volume Vco{t) is modelled using Eq. (16.101) . A rather simplified 
prescription is adopted to treat the merging of correlation volumes: If at any given 
moment r the total correlation volume (sum over all "bubbles") kVco{T) is larger 
than Vfb{t), kVcoij) is set to Vfb{j) for the subsequent evolution. In order to check 
the sensitivity of J /ip regeneration to different sizes of the correlation volume, we 
multiply Vco{t) from Eq. (I6.10p with different overall scaling factors and compare 
the resulting inclusive yield. The results are presented in Fig. [541 First we see 
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Fig. 54. J /ip regeneration with different charm-quark correlation volumes. Left panel: 
thermal charm-quark spectra; right panel: longitudinal thermal + transverse 
pQCD charm-quark spectra with angular correlation between c and c. Dou- 
ble-dot-dashed line: regeneration in grand-canonical ensemble; dot-dashed 
line: regeneration in canonical ensemble, but without correlation volume ef- 
fect; solid line: regeneration with Ix V^o; dashed line: regeneration with 2xVco'i 
dotted line: regeneration from the thermal rate equation approach. 



that the charmonium regeneration is stronger in the canonical ensemble than in the 
grand-canonical ensemble. Including the correlation volume effect further enhances 
the regeneration due to the effectively larger probability for one charm quark to find 
its partner (by a factor of VpElkVco)- Because Vco grows faster than Vfb this effect 
is more pronounced for peripheral collisions where the medium lifetime is relatively 
shorter and the early stage regeneration has greater impact on the final yield. We 
also note that with Vco doubled to 2Vco the regeneration yield almost drops by a 
factor of two (dashed line). However, further increasing Vco the regeneration yield 
only decreases by a limited amount due to the fact that for large correlation volume 
kVco ~ Vfb the correlation volume "bubbles" begin merging and their maximal size 
is restricted to the entire fireball volume, Vfb- This can be seen from the fact that for 
central collisions the regeneration with 2Vco is already quite close to the limit (Dot- 
dashed line) where the cc pairs are correlated inside the entire fireball (essentially no 
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correlation volume effect). With the correlation volume effect included the inclusive 
yields of J /if) regeneration from thermal and from longitudinal thermal + transverse 
pQCD charm-quark spectra (with cc angular correlation) span a rather large range. 
The rate equation result (with 1 x Vco) is inside these two limits. 

Our studies in this section demonstrate the fact that J /ip regeneration is very 
sensitive to the charm quark off-equilibrium effects in both momentum (partial ther- 
malization) and coordinate space (correlation volume effect). It is therefore of crucial 
importance to implement these effects in a more systematic and more realistic way. 
Ideally one would obtain the time-dependent joint c and c phase space distribution 
from, e.g., Langevin simulations, as input to address these issues. Work in this direc- 
tion is planned. 

D. i})' and Xc Production 

In addition to providing a feeddown contribution to J /ip production, excited charmo- 
nia can give valuable complementary information on the medium created in heavy-ion 
collisions. Being more loosely bound states they usually have larger dissociation rates 
than J /if), and due to their heavier masses, their thermal equilibrium abundances are 
smaller than J /ip. Therefore we expect stronger suppression for excited charmonia 
compared to J /ip. However, if formation time effects |128[ I130[ 1129] are important, 
one may observe less suppression for Xc than for J /ip. 

With the thermal rate equation approach we calculate the ip' to J /ip ratio in 
y^=17.3A Pb+Pb collisions and compare to NA50 measurements in Fig. [551 In 
p+p collisions at SPS energies, the ratio of produced ip' to J /ip mesons amounts 
to a value of about 0.017 (the branching ratios from ip' and J /ip into dimuons are 
included). The nuclear absorption cross section for ip' , a^j^g, extracted by NA50 is 
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Fig. 55. Inclusive ip' to J/ip ratio at SPS compared to NA50 data |137l 1138] . Solid 
line: strong-binding scenario; dashed line weak-binding scenario. 

7.9 mb compared to aJ^f=4:Ainb for J/4'- Both these two numbers are measured in 
400 GeV p+A collisions. With the updated cr^,^g =7.3 mb at 158 GeV and assuming 
that the ratio between af-^^s and a^lJ is the same for p+A collisions at both 158 GeV 
and 400 GeV, we obtain a^^^3=13mb for v^=17.3AGeV Pb+Pb colhsions. The data 
from NA50 suggest additional suppression in ip' relative to J/ip in Pb+Pb collisions 
on top of the larger nuclear absorption cross section for ip'. In the strong-binding 
scenario this can be explained by the larger quasifree dissociation rate in QGP due 
to the smaller binding for ip'. In the weak binding scenario a significant fraction of 
charmonia is from regeneration, less ip' than J/ip are regenerated due to the lower 
dissociation temperature of ip', T^^^ ~ Tc vs. TJ^^^ ~ 1.25Tc, implying a later onset 
of regeneration for ip'. 

For \1/ production in A/i=200AGeV Au+Au collisions at mid-rapidity at RHIC 
we use a^l^=a^^^=3.5mb and aj^^^ ~6.5mb. The ip' to J/ip ratio and the x to J/ip 
ratio are displayed in Fig. |56l For Xc states, we constrain ourselves to xa and Xc2 
with a combined average branching ratio of 27% into J/ip^s. Both drop with centrality 
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Fig. 56. Inclusive {xa + Xc2) / {J / i^) ratio (left panel) and ip' / {J / ip) ratio (right panel) 
vs. centrality at mid-rapidity at RHIC evaluated with the thermal rate-equa- 
tion approach. Solid line: strong binding scenario; dashed line: weak binding 
scenario. Dotted line: nuclear absorption only. The J/^'s in the denominator 
include feeddown from Xc and ijj' . 

below the ratios obtained from CNM-induced suppression. 
E. Charmonium Production at FAIR 



Finally we briefly discuss charmonium production at the forthcoming FAIR accelera- 
tor. The typical collision energy for heavy nuclei is up to ^i ~10 AGeV. The medium 
created at FAIR is expected to have lower initial temperature and larger baryon den- 
sity than at SPS, recall Figs. |30]and[3ll In central collisions QGP is still expected to 
form. 

We calculate charmonium production in y^=8.8AGeV Pb+Pb collisions using 
the thermal rate equation approach. For the charmonium production cross in p+p 
collisions we take d(T^/(iy=2.4nb [91]. We use 0"^^^^ =7.3mb |106j and cr!fjjg=13.0mb 
by assuming the nuclear absorption at FAIR is the same as at SPS. The input charm- 
quark cross section in p+p collisions is taken as dcrgc/dy (y=0)=0.2/ib according to 
the extrapolation in Ref. |139] . 
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Fig. 57. Temperature dependence of J /if) binding energy at FAIR in the strong binding 
scenario. It is obtained by rescaling the J /ip binding energy at /i_B=0 according 
to Eq. (EH]). 
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Fig. 58. Results of the thermal rate-equation approach for i?y^ vs. centrality at FAIR. 
Solid line: total J /ip yield; dashed line: suppressed primordial production; 
dot-dashed line: regeneration component; dotted line: primordial production 
with nuclear absorption only; double-dash-dotted line: suppressed primordial 
production with quasifree suppression in the QGP only (without hadronic 
suppression) . 
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Since currently IQCD can only yield reliable results in the low baryon density 
regime, we lack information on in-medium \E' properties from first principle calcu- 
lations for FAIR conditions. For exploratory studies we obtain the in-medium "^ 
binding energies at FAIR by rescaling those at RHIC energy according to the ratio 
between the critical temperatures, Tj^^^ and T^^^'^, namely. 



ei?(7')|FAIR - (^B ( ^FAIR ^ 



Z 



RHIC 



' c 



(6.14) 

RHIC 

where T^^^ic jg around 180 MeV and T^?^^^ is around 135 MeV, recall Fig. |3T1 In 
this dissertation we restrict ourselves to the strong binding scenario. The resulting 
in-medium binding energy for J/ip is illustrated in Fig. [571 

We present the resulting centrality dependence of J/ip Raa in Fig. [SHI For central 
collisions the regeneration component takes up a significant fraction and compensates 
for the anomalous suppression of the primordial component, rendering the total yield 
even slightly above that resulting solely from CNM-induced suppression. Most of the 
suppression of the primordial component is due to nuclear absorption. The majority 
of the anomalous suppression is from partonic dissociation in the QGP. The suppres- 
sion in hadronic matter is negligible. Note however that the SU{A) effective theory 
employed here for hadronic dissociation takes into account only mesons-induced sup- 
pression, while the medium at FAIR energy is rather baryon-dense. Therefore the 
hadronic suppression shown in Fig. [58] must be considered as a lower limit, susceptible 
to significant corrections from inelastic \E' collisions with baryons as studied, e.g., in 
Ref. [HO]. 

Next we evaluate the centrality dependence of the ip' to J /ip ratio. In p+p 
collisions in the FAIR energy regime this ratio is extrapolated to be around 0.01 [1^ 
(The branching ratios from ip' and J/ip into dimuons are included). The results of 
the rate-equation approach, presented in Fig. [591 show that the ■?/''/( J/ ■?/') drops with 
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Fig. 59. Results of the thermal rate-equation approach for the tp' /{J/tp) ratio vs. cen- 
trality at FAIR. Solid line: full result; dotted line: primordial component with 
CNM effects only. 
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Fig. 60. Results of the thermal rate-equation approach for (p^) vs. centrality at FAIR. 
Solid line: total J/^ yield, dashed line: suppressed primordial production; 
dot-dashed line: thermal regeneration; dotted line: primordial production 
with CNM effects only. 
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centrality well below the ratios obtained from CNM-induced suppression as a result 
of anomalous suppression. 

We finally turn to J /ip transverse momentum spectra. The {Pt)pp in P+P colli- 
sions at FAIR energy regime is extrapolated to be around 0.8 GeV^ |142j . We assume a 
Cronin effect at FAIR similar to that of SPS, with agAr=0.076GeV^/fm. The resulting 
centrality dependence of {p^) is displayed in Fig. (HHl Similar to SPS energies the (p^) 
for the primordial and regeneration component exhibits a different centrality depen- 
dence, where the former is largely determined by the Cronin effect. The anomalous 
suppression in QGP induces a small suppression of {pf) due to the larger dissociation 
rate for \E' with higher momentum, recall the right panel of Fig. [211 The (pf) for the 
regeneration component increases with centrality due to the growing collective flow. 
The pt dependence of RJ^ for selected centralities is summarized in Fig. [6TJ Again, 
the Cronin effect prevails via the primordial component and the collective flow entails 
the Raa of the regeneration component to increase with pf. 



131 



1.5- 



J^ 1- 



0.5 



■ Pb-Pb (8.77 A GeV) N =365 

a^j^=0.076 GeVVfm 


■ 


- ■ ■ ■ Nuc. Abs. 

— total 

— primordial 


/: 


■ - regeneration 


/ 


___^ 







1.5 



- J^ 1 



0.5 



Pb-Pb (8.77 A GeV) 

a.„=0.076 GeV7fm 

Nuc. Abs. 
total 

primordial 
regeneration 




0.5 1 1.5 

P, (GeV) 



0.5 1 1.5 

P, (GeV) 



1.5 



< 1 - 



0.5 



Pb-Pb (8.77 A GeV) 


N =72 

part 


- a,_^,=0.076 GeV/fm 


, 


- ■ ■ ■ Nuc. Abs. 


/ 


— total 


/ 


— primordial 


/ 


■ - regeneration 


y^ : 


^. ^ . ^ ., 1. ^ . ^ . ^ ., 


- 



1.5 



< 1 



0.5 



0.5 1 1.5 

P, (GeV) 



Pb-Pb (8.77 A GeV) 


N =25 

part 


a „=0.076 GeV7fm 

gN 


- 


- ■ ■ ■ Nuc. Abs. 


- 


— total 


/ 


primordial 


/ 


■ - regeneration 


^ ; 




- 


" 


1 



0.5 1 1.5 

P, (GeV) 



,J/t^ 



Fig. 61. -R^^ vs. transverse momentum for different centrality selections of Pb+Pb 
collisions at FAIR. Solid line: total J jtl) yield; dashed line: suppressed pri- 
mordial production; dot-dashed line: thermal regeneration; dotted line: pri- 
mordial production with CNM effects only. 



132 



CHAPTER VII 

CONCLUSIONS AND OUTLOOK 
Quantum Chromodynamics (QCD) has been established as the underlying theory for 
the strong interaction decades ago, but several open problems persist in the non- 
perturbative regime of the strong interaction to date. One of the prominent problems 
is to understand the phase structure of the QCD matter which ultimately consists of 
quarks and gluons. At low temperatures quarks and gluons are confined into hadrons, 
which become the effective degrees of freedom of the QCD matter. First principle 
numerical calculations of discretized (lattice) QCD predict that at a temperature of 
~170MeV (lO^'^K), ordinary hadronic matter will undergo a transition into a decon- 
fined phase where quarks and gluons are the relevant degrees of freedom, forming the 
so-called Quark-Gluon Plasma (QGP). One expects to gain insights into the strong 
interaction and many-body dynamics of QCD through investigating the properties of 
QGP. However, since QGP exists only at very high temperature, its natural occur- 
rence is rare. Fortunately QGP can be created in present-day's laboratories through 
ultra-relativistic heavy- ion collisions (URHICs). Since in URHICs QGP exists for a 
very short time (several fm/c), its properties can only be inferred from finally ob- 
served particles. Among these particles the charmonia as the bound states of charm 
(c) and anti-charm (c) quarks turn out to be an excellent probe for QGP in URHICs: 
The strong force binding c and c in vacuum is expected to be screened in QGP by 
surrounding colored quarks and gluons. As a result, charmonium bound states are 
easier to dissociate by collisions with particles in the medium, leading to a reduction 
of its experimentally observed yield. 

In order to utilize charmonium as a quantitative probe for the properties of the 
medium created in URHICs, one needs a framework associating experimental observ- 
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ables with theoretical calculations of in-medium charmonium properties. In this work 
we have constructed such a framework in which a Boltzmann transport equation is 
employed to describe the time evolution of the charmonium phase space distribution in 
URHICs. The in-medium charmonium properties figuring into the transport equation 
are constrained by thermal lattice QCD (IQCD): We have estimated the charmonium 
in-medium binding energies from a potential model with input 2-body potentials in- 
ferred from the free energy of a static pair of heavy quark and antiquark computed in 
IQCD. Based on the obtained charmonium in- medium binding energies a "quasifree" 
approximation has been employed to evaluate the charmonium dissociation and re- 
generation rates, which are the main inputs to the transport equation. Finally, the 
consistency between the charmonium binding energy and dissociation rates has been 
verified with the quarkonium current-current correlation function, another quantity 
calculated by IQCD independently from the heavy quark free energy. In this way we 
have established a link between the equilibrium properties of charmonia calculated 
from first principles, but in euclidean spacetime, and the off-equilibrium evolution of 
charmonium phase space distributions in heavy-ion collisions. 

In the Boltzmann equation the dissociation and regeneration rates play the role of 
the interface between the microscopic dynamics and the macroscopic observables. In 
this dissertation we have worked out the momentum dependence of the charmonium 
dissociation rate, which turns out to be sensitive to the charmonium binding energies. 
We have also explicitly worked out the charmonium regeneration rates based on the 
detailed balance between the dissociation and regeneration processes, not restricted 
to equilibrium conditions. We have found that the dominant regeneration process in 
QGP is a 3-to-2 process with the initial states consisting of c, c and a light parton. 
In the "quasifree" approximation this 3-to-2 process is found to be factorizable into 
a 2-to-2 scattering process and a 2-to-l coalescence process. We have calculated the 
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regeneration rates with different input charm quark momentum spectra. We have 
found that, the thermal charm-quark spectra turn out to be much more efficient in 
regenerating charmonia than the charm-quark spectra from initial hard collisions, due 
to a larger phase space overlap between c and c quarks. Both the inclusive yield of 
regenerated charmonia and their pt spectra appear to be very sensitive to the level of 
thermalization of the input charm quark spectra. 

Within the current uncertainties from various inputs our results from the trans- 
port equation agree reasonably well with the J/ip production data measured at the 
Super-Proton-Synchrotron (SPS) and the Relativistic Heavy-Ion Collider (RHIC), 
thus corroborating the picture of the deconfining phase transition as predicted by 
QCD. Moreover we have demonstrated that the J/ip transverse momentum data 
exhibit promise for discriminating power for different scenarios for the heavy-quark 
potential as extracted from IQCD data: The scenario with the internal energy (rather 
than the free energy) of the cc system identified as the potential is slightly favored, 
though further theoretical studies and more precise experimental data are required 
to draw definitive conclusions. 

To reduce the current uncertainties of the theoretical approach, further devel- 
opments are in order. First, our calculations show that in URHICs the partition 
between the primordial and regenerated charmonia is sensitive to their in-medium 
binding energy, therefore a reliable determination of the in-medium c-c 2-body po- 
tential is needed for an accurate evaluation of charmonium in-medium binding energy, 
and thus of the dissociation and regeneration rates. 

Second, an explicit calculation of charmonium regeneration from (time-dependent) 
charm-quark phase-space distributions in c-c recombination reactions should be per- 
formed. One should use realistic charm-quark spectra as following, e.g., from Langevin 
simulations with constraints from the T-matrix formalism and from open-charm ob- 
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servables, to reduce the currently large uncertainty in this part of the input. 

Third, hydrodynamic simulation of the medium evolution could be employed for 
a more detailed and realistic description of the temperature and the flow field of the 
underlying medium, especially in coordinate space. 

Fourth, a microscopic model for primordial cc and charmonium production is war- 
ranted to better disentangle nuclear shadowing and absorption in the pre-equilibrium 
stage, including formation-time effects. This would improve the initial conditions for 
the transport approach in the hot medium. 

These developments will ultimately lead to a comprehensive approach which 
can serve as a quantitative bridge between charmonium phenomenology in heavy-ion 
collisions and theoretical studies of charmonia in the QGP (and hadronic matter). It 
will enable us to deduce from experimental observables effective degrees of freedom 
for the hot and dense medium created at various collision energies, improving our 
knowledge of the phase structure of hot and dense QCD matter. Furthermore it will 
provide insights into basic properties of the strong force in terms of color screening of 
Coulomb and confining interactions. All these aspects of information will eventually 
contribute to the establishment of a coherent picture of the strong interaction in the 
non-perturbative regime. 
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